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ABSTRACT 


The theory of acoustic plane wave normally incident on a clamped 
panel in a rectangular duct is developed from basic theoretical concepts. 
The coupling theory between the elastic vibrations of the panel (plate) 
and the acoustic wave propagation in infinite space and in the rec- 
tangular duct is considered in detail. The partial differential equa- 
tion which governs the vibration of the panel (plate) is modified by 
Adding to it stiffness (spring) forces and damping forces, and the 
fundamental resonance frequency f and the attenuation factor a are 
discussed in detail. 

The noise reduction expression based on the present theory is 
found to agree well with the corresponding experimental data of a 
sample aluminum panel in the mass controlled region (f > f ) , the damping 
controlled region (f ~ f Q ) and the stiffness controlled region (f < f Q ) . 
All the frequency positions of the upward and downward resonance spikes 
in the sample experimental data are identified theoretically as resulting 
from four cross interacting major resonance phenomena: the cavity 

resonance, the acoustic resonance, the plate resonance, and the wooden 
back panel resonance. Detailed tables are given for the values of 
these resonance frequencies in each case. 
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CHAPTER I 


INTRODUCTION 

The acoustic research project started on April 15, 1976, when the 
Flight Research Laboratory of the University of Kansas b;egan work under 
a grant from NASA, Langley Research Center, entitled, "A Research Program 
tjo Reduce Interior Noise in General Aviation Airplanes" (NASA Grant 
No. 1301) . Over the past three years a research facility has been 
established, in the form of a Beranek tube and additional equipment, 
and a large volume of experimental data has been published. A sub- 
stantial amount of the experimental data published to date involves 
the noise reduction curves for acoustic plane waves normally incident 
on clamped aluminum panels, and other materials, in the Beranek tube. 

The aim of the present report is to present a theory which will 
explain reasonably well the detailed features of the experimental noise 
reduction curves for normal incidence of acoustic plane waves on a 
clamped panel in a rectangular duct (the Beranek tube) . Such detailed 
features include the general behavior of the noise reduction curve in 
all its parts, as well as the frequencies of the numerous resonance 
spikes, which are superimposed on the curve both upward and downward, 
for frequencies above the fundamental frequency of resonance. Another 
aim of the present report is to be used as a theoretical base on which 
more complicated acoustics experimental results could be studied from a 
theoretical point of view in the future. 

Such a theory sould be based on the interaction between two 
general fields of study: The theory of acoustic wave propagation in 
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infinite space and in ducts, and the dynamic theory of plates and the 
theory of elasticity. Since the vibrations in the panel (plate) and 
the acoustic waves in the air are coupled strongly along the whole 
panel, and are affecting each other noticeably, the interaction between 
these two systems plays a major role in this theory. When one sends an 
incident wave of sound in the direction of the panel, this wave will 
be reflected from the panel. At the same time it will set the panel into 
motion of vibrations, which will generate transmitted acoustic waves 
on the other side of the panel. The incident, reflected and transmitted 
acoustic waves will be coupled strongly together with the induced 
transverse displacement of the panel. Because of the experimental 
set up in the Beranek tube, the acoustic waves propagate in a duct with 
a square cross-section. Higher order acoustic modes, as well as the 
fundamental mode of plane waves, should be taken into account, in 
order to explain the experimental results of resonance at certain 
particular frequencies . 

The present report develops the theoretical derivations from the 
basic equations. The theory presented in this report successfully 
explains the detailed features of the noise reduction curves for one 
particular aluminum panel, which has-been taken as a typical example. 

It should be pointed out that one of the most challenging aspects of 
the present theory, which has been met successfully, is to avoid lengthy 
numerical analysis, in order not to mask the main features of the 
interaction between the different physical phenomena. All the calculations 
in the present report have been done on a simple hand calculator. 
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The Meter, Kilogram, Second system of units is being used through- 
out this report , except at some places where the experimental data is 
given otherwise. The factor e la)t is used for harmonic time variation. 

In this report the following vector identities will be used, where 
* x » a y , a 2 are the unit vectors in rectangular coordinates: 


Vp = 

V • 

V x u 


(a — + a -$-+ l ~) p = a -^ + I + 2 

x 3x y 3y z dz ^ x 3x y 9y z 3z 


a 


; 4- + s 

-L + i J_) 

3y a z 3z' 

• u 

3u 3u 3u 

* + 

x 3x 

y 


3x 3y 3z 

3u 

3u 

3u 

3u 

3u 3u 

~r~ - 


— ) 
3x ' 

+ s z ( lf - 

dy 

3z 

y 3z 


V 2 p = 7 • 7p = + A 

3x 3y 3z 


V x Vp - 0 V • Vx u = 0 

— - 2 - 

V x (V x u) = V(V • .u) - V u 



closed surface integral 


closed surface integral 
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CHAPTER II 


ACOUSTIC WAVES - BASIC THEORY 

In the present chapter the basic equations of the acoustic waves 
will be given to be used in this report. 

The equation of state for a perfect gas may be given in the form: 


MP - R pT (1) 

o 

where: M = molecular weight of the gas ( Kg/mole) 

2 3 

P = the total pressure of the gas (N/m ) or (J/m ) 

3 

R = 8.314 x 10 ( J / °Kmole)-gas constant 

o 

3 

p = the total gas density of the gas (Kg/m ) 

T = temperature in °K where T(°K) - 273.15° + T(°C) 


In the particular case of air, which consists of Oxygen (M 
and Nitrogen (M^ = 28 Kg/mole), the molecular weight is: 


o 2 


= 32 Kg/mole) 


M . = 28.9644 (Kg/mole) 

3.X IT 

Substituting this value in Eqn. (1) one obtains the equation of state 
for air in the form: 


P = RpT (2) 

where the gas constant for air is given by: 

R 

R = = 287.05 (J/°K Kg) 
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For example, the U.S. standard atmosphere tables give the following 
values at sea level at T = 15°C * 59°F: 

P = 760mm Hg = 1.01325 x 10 3 (N/m^) (1 atmosphere = 10 3 N/m 3 ) 
p = 1.2250 (Kg/m 3 ) 

T= 273.15° + 15.00° = 288.15° 


and obey Equation (2) above. 

The acoustic wave motion, when the wave amplitude is small, can 
be described by the following varying quantities: 


2 

p =■ p(x,y,z,t) = acoustic or sound pressure (N/m ) where |p| << P 
u - u(x,y,z,t) = velocity vector of the air (m/sec) 


where P is the equilibrium pressure. 

One requires one scalar equation and one vector equation to relate these 
varying quantities for the acoustic wave motion. 

From the equation of continuity one obtains: 


K 3p(x,y,z,t) 


= -V • u(x,y, z, t) 



(3) 


which states that the velocity gradient produces a compression of the 
air. The only forces involved are those of compressive elasticity, 
measured by the compressibility k. Since in the case of the acoustic 
wave motion there is zero heat conduction in the air, one should use the 
adiabatic compressibility given as follows: 
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1 2 3 

k = yp = adiabatic compressibility (m /N) or (m /J) 

where P is the equilibrium pressure and for perfect diatomic gas, such 
as air, one has; 

C 

s 2 * 1 4 s _£ s specific heat at constant pressure 
Y 5 C v specific heat at constant volume 

From the equation of wave motion one obtains; 


o 8 u(x,y.z,t) 
St 


vp(x,y,z,t) 


(4) 


which states that a pressure gradient produces an acceleration of the 
air, where p is the equilibrium density of the air. 

Taking the divergence on both sides of Eqn. (4) and substituting 
Eqn. (3) one obtains: 

2 

V • Vp = (v • u) = -p-^(-ic|£) = <P ~ 2 

3 1 

which may be rewritten in the form of the wave equation: 

2 2 

V 2 p - icp— & = 0 or V 2 p — t = 0 (5) 

st c at 2 

where c is the velocity of the acoustic wave in the air. Substituting 
K = ^ in Eqn. (5), using Eqn. (2), one obtains: 

— .=.A£*/ yRT = 20.047/T(°K) 

/kp V p 


c 


(m/ sec) 


( 6 ) 


where one has: 


T(°K) * 273.15° + T(C°) 

For example at T * 15°C (59°F) the acoustic wave velocity is: 

c = 340.3 m/sec for T = 15°C(59°F) 

Taking the gradient on both sides of Eqn. (3) and substituting Eqn. (4) 
one obtains: 

V7 • G = -<c^(Vp) = [-pf£] = 

9 1 

Since from Eqn. (4) one has for the time varying acoustic wave V x u = 0, 
one obtains: / 

2- 

VV • u - V u x V x u = 0 since V x u = 0 

v 

Substituting it in the above, one has the wave equation for the veolcity 
vector field: 

2- 3^u n „2- 1 8 2 u n ,_ s 

V a - Kp — ~ =0 or V u r ~ = 0 (7) 

3t c 3t 

In order to find the transmitted power in the acoustic wave motion 
it is necessary to develop a power theorem for the acoustic waves, 
using the vector identity: 
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V • (pu) « pV i u + u • Vp 


( 8 ) 


By multiplying both sides of Eqn. (3) by p and dot-product both sides of 
Eqn. (4) by u and substituting in Eqn. (8) one obtains: 


V • (pu) = -kp|E- - pu • 


(9) 


It may be shown that! 


9P 13/ \ 9/i 

Kp 9t *2 IF (pp) = > 


3 /I 2 n 


( 10 ) 


- 3u 1 3 -v 3 /1 | - 1 2v 

PU * 9l = p 2 9T (u * U) = H ( 2 p|u l ) 


9 ,1 2 % 

9l ( 2 pu } 


(ID 


Substituting Eqn. (10 and Eqn. (11) in Eqn. (9), one obtains: 


V • (pu) = - ^(|kp 2 ) - ^(|pu 2 ) (12) 


Integrating both sides in an arbitrary volume V and using the Gauss 
divergence theorem, one obtains 



V • (pu)dV 



ndS = - 



1 2 

+ ypu]dV 


(13) 


Physically Eqn. (13) states that the total radiated (transmitted) 
acoustic wave power through the surface S is equal to the total decrease 
of acoustic wave power inside the volume V, where we have: 
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1 2 
w p * 2 Kp 

1 2 

w R = pu 
I = pu = 


3 

■ acoustic wave potential energy density (J/m ) 

3 

= acoustic wave kinetic energy density (J/m ) 

instantaneous radiated (transmitted) acoustic wave 

2 2 
power density (Watt/m ) or (J/sec m ) 


Thus the radiated (transmitted) acoustic wave power is in the direction 
of the acoustic wave velocity vector u. The vector I is sometimes 
called the "sound intensity." Taking the total acoustic wave energy 
to be w = Wp + Wj,, one can rewrite Eqn. (12) in the form: 


V . i + “(Wp + w R ) = v’ • I + ~w = 0 (14) 

where Eqn. (14) gives us the nondissipative energy equation of continuity 
for the acoustic wave with no sources in the vicinity. 

By using the Gauss divergence theorem, Eqn. (3) may be rewritten 
in an integrated form as follows : 



dt 


ffh ‘Iff ' “ dV # 


ndS 


V 


(15) 


where the right hand integral gives the total velocity flux. 

By taking a small cylinder of area (AS) and height (Az) across the 
boundary between two different fluid media, applying the integrals in 
Eqn. (15) and taking Az -*■ O.one obtains the following boundary condi- 
tion: 



( 16 ) 
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The normal component of the velocity vector of the acoustic wave is 
continuous across the boundary. 

By using the corresponding vector theorem, Eqn. (4) may be 
rewritten in an integral form as follows: 



# 4V ■ Ilf 


Vp dV = 


(17) 


where the! integral /X/pudV gives the total momentum in the volume. 

V 

Similarly to the above, one obtains the following boundary condition: 



(18) 


The pressure of the acoustic wave is continuous across the boundary. 
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CHAPTER III 


ACOUSTIC UNIFORM PLANE WAVE 

In the present chapter the acoustic waves basic equations will be 
applied to the harmonically varying acoustic uniform plane wave in 
infinite fluid (air) medium where the effects of viscosity, heat con- 
duction, relaxation, vibration and other dissipative (attenuation) 
processes have been neglected. 

Assuming harmonic time variation e la)t , where m is the circular 
frequency, one may write the basic quantities of the acoustic wave by 
using phasor convention as follows: 


p(x,y,z,t) = Re/2p(x,y,z)e lu)t 

(19a) 

u(x,y,z,t) = Re/2 u(x,y,z)e iu>t 

(19b) 


where Re represents the real part and p(x,y,z) and u(x,y,z) are the 
complex effective value phasors, which represent both amplitude and 
phase in the complex plane. For example one has: 



p = Re/2 1 p | e^e =/2 |p| cos (cot - 4>) 


where |p| is the effective value, with similar relations for each 
component of the complex velocity vector u(x,y,z). 

In order to simplify the notation, p (x, y, z) and u(x,y,z), will be denoted 
from now on as p = p(x,y,z) and u = u(x,y,z), where the harmonic time 
variation convention of Eqn. (19) is understood. 
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1(jO t 3 

Assuming harmonic time variation e and replacing — - by (-iw) 
the two basic equations Eqns . (3) and (4) for the acoustic waves for 
harmonic time variation will becomes 


> 


iojKp = V • u 


(20a) 


iwpu = Vp (20b) 

where complex quantities p = p and u = u in Eqn. (20) are understood. 

From (20b) one also obtains: 

V x u = 0 (20c) 

Thus, u(x,y,z) may be found from a given p(x,y,z)in Eqn. (20b) and 
p(x,y,z) may be found from a given u(x,y,z) in Eqn. (20a). By taking the 
divergence of both sides of Eqn. (20b) and substituting Eqn. (20a), one 
obtains : 

2 

V • Vp = icopV • u = iu)p(icotcp) = — (uj Kp)p 

which may be written in the form of the wave equation for harmonic 
time variation • 

V^p + (m^Kp)p = 0 or V 2 p + k 2 p = 0 (21a) 
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By taking the gradient of both sides of Eqn. (20a) and substituting Eqn. 
(20b) one obtains 

2 - 

V(V • u) = iuiicVp = iu>ic(iiDpu) = -(u> kp)u 

By using Eqn. (20c) it may be written in the form of the wave equation 
for harmonic time variation: 

V 2 u + (u) 2 Kp)u = 0 or V u + k u = 0 (21b) 

The wave number k is defined by: 

k = a)i^Kp"= ~ (l/ m ) (22) 

where: w = circular frequency (l/sec) 

7 T- = f = the frequency of the acoustic wave (l/sec) 

2ir 

X = the wavelength of the acoustic wave (m) 

f\ - c = the velocity of the acoustic wave (m/sec) 


Assuming acoustic uniform plane wave propagating in the z-direction 
one has p = p(z) and Eqn. (21a) becomes: 


dz 2 


+ k p = 0 


P = P Q e 


ikz 


or 


p = p Q e 


-ikz 


where p is a constant, 
o 

Substituting the first solution in Eqn. (19a) one obtains: 
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(23a) 


p(x,y,z, t) = Re/2 p e *(k z- wt) a/2 p. cos(kz-ut) 

o o 


Substituting Eqn. (23a) in Eqn. (20b), one finds u = 0, u ■ 0, 

x y 

u ^ 0. By solving Eqn. (21b), one obtains similarly: 
z 


u (x,y,z,t) = Re/2u =/2 u cos(kz-u)t) 


(23b) 


2_lr 2 lie Z 

Substituting p = p Q e and ~ u Q e in Eqn. (20a) and Eqn. (20b) 
one obtains by using Eqn. (22) : 


p _ j_o _k_ up. _ /p _1_ _ pc _ 


U U U)K 

z o 


(23c) 


where Z is the characteristic acoustic wave impedance in infinite 
o 

3 2 

medium in (Nsec/m ) or (kg/m sec). 

The characteristic acoustic wave impedance for air at sea level at 
15°C(59°F) will be: 


Q Wcpn 

Z = pc = (1.2250 Kg/ni ) • (340.3 m/sec) = 417 
o o 

m 


Equations (23) represent the acoustic uniform plane wave propagating 
in the positive z-direction. The pressure is in phase with the fluid 
velocity. It is a longitudinal wave since the velocity vector u has only 
one component in the direction of propagation. 

Similarly, the acoustic uniform plane wave propagating in the 


negative z-direction will have to form: 
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p(x,y,z,t) = Re/2 p Q e i ^ z+ w t )_ /2p o cos(kzfa)t) (24a) 

u (x,y,z,t) * -Re/2u e *(kz+ojt)_ -/Ju cos(kz-Hi)t) (24b) 
z J o o 


P Q P c 

where u = — = — 
o Z pc 
o 


The velocity vector u is again directed in the direction of propagation 


of the acoustic plane wave. 

The average acoustic plane wave potential energy density is 
given by: 


x 

”p = -b / [/?| 

t=0 


[/2p (? cos(kz-a)t)] 2 dt = p Q “ - f | p 1 2 (J/m 3 ) (25a) 


2tt 

where the integral has been evaluated by taking oj = 277 f - — 
where T = the period (sec) 

Similarly , the average acoustic plane wave kinetic energy density is 
given by: 


_ p 

W K 7T 


i 

j [ /2u o 

t=0 


cos'(kz-ojt)' -}-dt = "2 u ^ = Y l u (25b) 


2 2 

From Eqn. (23c) we have kP q = pu Q and therefore in the acoustic 

. . . . 

uniform plane wave one has w = w • The average total acoustic wave 

r K 

energy density is given by: 


w = w p + w K = k|p| 2 = p |u z | 2 = 

pc 


(J/m 3 ) 


(25c) 
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The average radiated (transmitted) acoustic plane wave power (the 

"sound intensity") is given by: 

• T 

X z = 2p o U o T /cos 2 (kz- u t)dt = p o u o = M- = pcjuj 2 * 
t=0 

= cw = -hp u + pu *) (Watt/m 2 ) (26) 

z z z 

* 

where p is the complex conjugate of p. 

Let an acoustic plane wave propagating in the positive z-direction 
be incident on a rigid wall at z = L. As a result an acoustic wave will 
be reflected from the wall propagating in the negative z-direction. 

The incident plane wave is given by using the phasor convention as 
follows: 



Px e 


ikz 


u 


zi 


1 xkz 
— D_e 
pc ‘ I 


(27a) 


where p,j. is a complex constant. 

The reflected acoustic plane wave is given by: 



-ikz 




1 -ikz 


(27b) 


where p^ is a complex constant. 

K 

At the rigid wall z = L one has u^ = 0: 



r . i 1 r ikL -ikL. n 

[u . + u ] = — [p T e - Pr, e ] = 0 

L zi zr T pc r I r R 

z=L 


(28a) 
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which gives: 


P R " P i e 


i2kL 


(28b) 


The total pressure and velocity fields of both waves will be therefore: 


, ikz . i2kL -ikz 

p = p i + p r = Pj e + Px e e 


(29a) 


1 ikz 1 i2kL -ikz 
u - u . + u = ■ — p_e - • — p T e e 

s zi zr pc I pc I 


(29b) 


where at z = L one has u = u = 0. 

z n 

-iojt 


Multiplying by /2 e <u and taking the real part of it in accordance 
withEqns. (19), one finds the actual physical fields: 


/ x , ikz . i2kL -ikz, -iu)t 

p(x,y,z,t) = Re/2 p^. [e + e e ]e 


- Re/2 p T [e i(Kz - kL) +e' i(kz - kL) J - 


= 2/2 PjCOs(kz-kL)cos(a>t-kL) 


(30a) 


. . _ 7x 1 r ikz i2kL -ikz, -icot 

u z (x,y,z,t) = Re/2 — p-^e - e e ] e 


- Re/2 -i- p T [e l(kZ ' kL) - e 1<kz - kL) ] 
pc I 


Re/2 — P x 2i sin(kz-kL) e i ^ u,t kL ^ = 
pc I 


2/2 — s in (kz-kL ) s in (ui t-kL ) 
pc 
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(30b) 


i *.r£ 


Both Eqn. (30a) and Eqn. (30b) represent standing waves. One may 
find the maximum value of the pressure from Eqn. (30a) to be at any 
point z = z . 


p (z=z ) = 2/2 p.lcos k(z -L) 
r max o r I' • o 


(31a) 


Eqn. (31a) may be found alternatively directly from Eqn. (29a) as 
follows: 


„(z=z ) = /2 


ikz n , i2kL -ikz ^ 1 
e 0 + e e 0 = 


= Jl Pjje 


ikL | | e ik(z °- L) + e“ ik(z o" L) |=2/2 p^cos k(z Q -L) 


(31b) 


Let us consider the case when the wall at z = L is not rigid. We 
will define the non-dimensional specific acoustic impedance £ of the 
wall at z = L as follows: 


pcu 


pc 


= ? = e - ix 


( d imens ionles s ) 


(32a) 


z=L 


Where 0 is the specific acoustic resistance and x is the specific 
acoustic reactance. Substituting Eqn. (27) in Eqn. (32a), one obtains: 


ikL , -ikL 
P I S + P R e 

ikL -ikL 

P l e “ P R e 


= ?; 


(32b) 


From Eqn. (32b) one obtains the complex reflection coefficient Cr: 


r P R 5-1 i2kL 
C r " ’ C+T e 


(32c) 


18 


For the particular case of a rigid wall ? = -i“ one obtains Eqn. 

For the particular case of a purely reactive surface 5 = -ix one 
has | C r | =1 and : 

r = _ -iX-1 i2kL _ -i2<|> i2kL _ i2(kL-<j>) 

r P r -ix+l 

The reflected wave differs only in the phase angle from the incident 
wave; the reflected acoustic power is equal to the incident acoustic 
power; therefore no energy is absorbed by the surface and one has 
standing waves as before. For the particular case of the rigid wall 
X = oo and (ji = 0. The total pressure and velocity fields of both 
waves will be therefore: 


(28b). 


(33) 


e " Pi + P r * 


ikz , i2(kL-<f>) -ikz 

r e + P T e e 


U 

Z 


u . + u 
zx zr 


1 xkz 1 x2(kL-<j>) -ikz 

— p e pe T e 

pc F I pc F X" 


One may find the maximum value of the pressure from Eqn. (34a) to 
be at any point z = z q . 


(34a) 


(34b) 


s \ 7 n i ikz Q . i2(kL-di) -ikz ol 

<z= 2 o> ■ /2 p ll e + * e °|- 

= /2 P |e i(kL-<f)) l | e i [ kz o-( kL -<|>)] + e -i[kz 0 -(kL-(j,)] j = 

= 2/2 p I |costk(z Q - L) + <()] | = 

= 2/2 p I |cos[k(L - z q ) - <j> ] | (35) 
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If the acoustic impedance of the wall is a function of frequency, one 
has z = ?(w) and x = x(“)» and one will have in Eqn. (35) $ = $(<*>)• 


CHAPTER IV 


ACOUSTIC WAVES IN RECTANGULAR DUCTS 

In the present chapter the guided acoustic waves in a rectangular 
duct (tube) with rigid walls will be developed. By a guided acoustic 
wave it is meant that the direction of the acoustic energy flow must 
be primarily along the direction of the guiding system of the duct. 

The analysis will be limited here to the guiding system of the duct 
which is straight and has a uniform rectangular cross-section and 
rigid walls . 

The basic equations for the acoustic wave motion are given in 
Equations (3), (4) and (5) as follows: 

icf£ = -V • u (36a) 

p~£ = -Vp (36b) 

2 

V 2 p ~ - 0 where c - — — (36c) 

c 3t /kp 


The acoustic wave propagating along the uniform guiding system of the 
duct in the z-direction can be described in terms of a propagation 
factor e^^ Z where we have: 


p(x,y,z, t) = Re/2 p(x,y)e^ qz (N/m 2 ) (37a) 

u(x,y,z,t) = Re/l u(x,y)e i ^ Z (m/sec) (37b) 
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where p(x,y) is a complex scalar phasor and u(x,y) is a complex vector 
described by three scalar phasors u x (x,y), u^(x,y) and u z (x,y). 
Substituting Eqn. (37) in Eqn. (36) one obtains: 

9u 9u 

ioiKp = + iqu g (38a) 


itopU x “ lx 5 itupu y “ ia) P u z = i<lP 

(38b) 

+ ^4 + (k 2 - q 2 )p = 0 

03c) 

9x 9y 



where p = p(x,y), u x = u x (x,y), u y = u y (x,y), u z = u z (x,y) and 

k = » = (l/n) (38d) 

C /— C LA A 

vicp 

k being the infinite medium wave number. By solving the wave Eqn. (38c) 
one obtains p, and substituting in Eqn. (38b) one obtains u. 

Substituting in Eqn. (38c): 

2 2 2 

ki = k' - q (39a) 

c 

one obtains: 

2 2 

•2-1 + + k j> - 0 (39b) 

•x 3y 2 c 


Eqn. (39a) may be rewritten in the form: 



where q is the duct wave number, k is the infinite medium wave number 

and k^ is the cut off wave number. When k - k^ one has q = 0 and 

there is no propagation. It should be pointed out that q could be either 

positive for guided waves propagating in the positive z-direction, or 

negative for guided waves propagating in the negative z-direction. The 

cut off wave number k is defined by: 

c 




f X 
c c 


(1/m) 


(40b) 


where f is the cut off frequency and X^ is the wavelength of an 
acoustic uniform plane wave at the cut off frequency f in infinite 
medium . One obtains from Eqn. (40a) : 

q = k /l-(^) 2 » k\f 1 - {-f) 2 = k\/l - (~ ) 2 (1/m) (40c) 

c 


since c = fX = f X . From Eqn. (40c) it is found that when f > f 
c c c 

the duct wave number q is real and the acoustic wave propagates without 
attenuation along the duct. When f < f^, the duet wave number q is 
purely imaginary and the acoustic wave attenuates exponentially and does 
not propagate in the duct. When f = f , one has q - 0 and there is 
neither phase shift nor attenuation along the duct and this condition 
is called the cut off of the mode. The duct effects the acoustic 
wave of a particular mode of propagation as a high pass filter. The 
acoustic wave of a particular mode propagates in the duct only if its 
frequency f is larger than a certain cut off frequency f such that 
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For the propagating range f > f £ the phase velocity Vp :b"f the 
acoustic wave mode propagating. :in the duct is given from Eqn. (38d) and 
Eqn. (40c) as follows: 


V = = 

p q 


C "" ~ — = — (m/sec) f > f £ (41a) 

\j i- i-f ) 1 \fi~- (-f-) 2 


where one has for the phase velocity V of the wave in the duct V > c. 

P P “ 

For the propagating range f > f the group velocity V of the acoustic 

c g 

wave mode propagating in the duct is given fron Eqn. (41a) as follows: 


V = m 


V 


g dq dV 

1 — — ■ * * 

V df 
P 


' ^ ( T > 2 ' V 1 • ( f )! 


(m/sec) f > ■ f c (41b) 


where V <. c. From Eqn. (41a) and Eqn. (41b) one has: 
S 


c = /V V 

P g 


(41c) 


The wavelength measured along the acoustic guided wave in the duct in 

the z-direction is the distance represented by a phase shift of 2 tt, and 

is denoted by X : 
g 


A = 

g q 


\J l-(-f) 2 \J 1 - (y-) 2 


(m) 


(41d) 


The wave number q for the propagating mode in the duct for f > f^ can 
be expressed as follows: 


2nf 


q ■ v" = V 


27rf 

fA 

g 


2tt 

A 

g 


(1/m) 


(42) 
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The acoustic wave impedance in the duct is defined similarly to the 
case of the impedance of the acoustic plane wave in the infinite medium 
Z Q from Eqns. (38b) and (40c) as follows: 


Z * -2- = — = — -- —22. = — ^ (Nsec/m^) or (Kg/m^sec) 

\Ti - <f > 2 \f ^¥ 2 

where _> Z q , For the propagating acoustic mode f > f , Z^ is a real 
number, and there is a transfer of average acoustic power in the duct. 
For the attenuating acoustic mode f < f , Z^ is an imaginary number, and 
there is no transfer of average acoustic power in the duct. 

Taking f^/f = cos0 one may obtain the following relations for 
the acoustic wave mode of propagation from Eqns. (40) , (41) , (42) and 
(43) : 

f c A 

-r- = ~ — = cos0 (44a) 

r A c 

q = k sin0 (44b) 


V = 


p sin0 


V_ = c sin0 
§ 


(44c) 




L 

sin© 


Z, = 


d sin0 


(44d) 


where one would consider the mode of propagation of the acoustic wave 
in the duct as an oblique multiple reflection of the acoustic plane wave 
between two parallel plane boundaries, where the angle of incidence is 
0 with the normal to the boundary. 


( 43 ) 
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The acoustic waves in the duct with the uniform rectangular 
cross section propagate in the z-direction. The boundaries (walls) 
of the inside rectangular cross section of the duct o.f dimensions "a" 
and "b" are defined by x = 0, x = a, y = 0 and y = b. Assuming that the 
boundaries (walls) of the duct are rigid, one should have the normal vel- 
ocity of the acoustic wave fluid to be u a 0 to each one of the walls. 

n 

Mathematically, the boundary conditions of the walls of the duct are 
expressed as follows: 

u = 0 x = 0; u = 0 x = a (45a) 

x x 



Using Eqn. (38b) in Eqn. (45), one may express the boundary conditions 
as follows : 


= 0 x = 0 

8x 

= 0 y = 0 


and x = a 


and y = b 


(46a) 

(46b) 


The partial differential Eqn. (36b) should be solved now, subject to 
the boundary conditions given in Eqn. (46). Using the method of 
separation of variables, one may assume a solution of Eqn. (39b) in the 
form: 


p(x,y) - X(x) Y(y) 


(47a) 
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Substituting it in Eqn. (39b) and dividing by X(x) Y(y), one obtains: 


_1_ dfxoo + _JL_ + k 2 = 

X <*> dx 2 Y ^> dy 2 c 


(47b) 


This equation is to hold for all values of x and y; since x and 
y may change independently of each other, then one should have each term 
to be a constant: 


1 d X(x) 2 

X , 2 x 

dx 


1 d 


hill - 


Y(Y) dy 2 


(47c) 


2 2 2 
k=k z + k 
c x y 


(47d) 


Solving Eqn. (47c), one obtains the general solutions: 


X(x) = A cos k x + B sin. k x 

X X 


Y(y) = C cos k^y + D cos k^y 


Applying the boundary conditions given in Eqn. (46) 5 one obtains the 
solutions in the form: 


X(x) '* A cos k x 
x 


k - — 
x " a 


(48a) 


Y(y) = C cos k^y 


k = — 
y b 


(48b) 


where m and n are independent positive integers or zero. 


From Eqns . (37a), (47a) and (48), one obtains the expression for the 


pressure field of the acoustic wave in the rectangular rigid duct with 


rectangular cross section: 
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(49a) 


, ■ v n m-rrx nirg i(qz-ajt) 

p(x,y,z,t) = Re/2 . p cos — cos —rf^ e ^ 

O ct D 


= /2 p cos cos cos (qz-cot) 

o a d 


By using Eqn. (49a) in Eqn. (38b), one obtains the expressions for the 
velocity vector field components of the acoustic wave in the rectangular 
rigid duct with the rectangular cross section: 


t . s „ pr . mir mux niry l(qz-ait) 

u (x,y,z,t) = Re/2 i — ~r sin —— cos — e M 
x ’ -- ojpa a b 


= -/2 -^L. sin cos sin(qz-ojt) 
oopa a b 


(49b) 


u (x.y.z.t) - Re^2 i -SI. cos SH. sin SJ2 e Kq^«t) 
y copb a b 


tx nm mirx . mry . , 

- -v2 — — cos sin — sm(qz-ojt) 

topb a b n 


(49c) 


u (x,y,z,t) = Re/2 p cos cos e^ pz 
z ■ top o a b 


_ rr P o mirx niry i(qz-tot) 

= Re/2 — - cos cos e ■ = 


/2 ~ cos cos .SIX cos(qz-tot) 

z d a b 


(49d) 


The Eqns. (49) give the pressure field and the velocity field of the 
acoustic wave of mode (m,n) in the uniform rectangular duct of rigid 
walls. The fields obey the boundary conditions in Eqn. (45) and Eqn. (46) 
The velocity field component u given in Eqn. (49d) is in phase with 
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the pressure field given in Eqn. (49a) and represents a transmission 
of average acoustic power in the z-direction for a propagating wave 
f > f in accordance with Eqn. (26) . The velocity field components 
u^ in Eqn. (49b) and u^ in Eqn. (49c) are 90° out of phase with the 
pressure field given in Eqn. (49a) and there is no transmission of average 
acoustic power in the x-direction or in the y-direction for a propagat- 
ing wave f > f since by Eqn. (26): 



J. 

/ 


cos (qz-ojt) sin(qz-ait) dt 


t=0 


cos2(qz-q)t) 

4Toj 


= 0 


t=0 


T 

J sin 2 (qz-ojt) dt = 
t=0 


since f = ^ and uT = 2irfT = 2-rr . 

For the case of f < f the mode (m,n) of the acoustic wave 
attenuates exponentially since q will be a purely imaginary number as 

seen from Eqn. (40c). In this case u in Eqn. (49d) will be also 90° 

; \ 

out of phase with the pressure p in Eqn. (49a) and there will be no 
transmission of average acoustic wave power in the z-direction. Thus 
in case f < f , the acoustic wave energy in this mode (ra,n) will not 
propagate, but will oscillate back and forth near the source of the 
acoustic wave power. 

From Eqns. (47d) and (48) one has for the acoustic wave mode 
(m,n) in the duct: 

k c ■ (k c>m,n -\/V + k / ' \/<f > 2 + <f )2 ^ < 50a ) 
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From Eqn. (40a) and Eqn. (50a) one has: 


q =^k 2 - k c 2 =^k 2 - (^) 2 - (^) 2 (^) 2 - (^) 2 (1/m) (50b) 


From Eqn. (40b) and Eqn. (50a), one has the cut off frequency f and 
the cut off wavelength for the acoustic wave mode (m,n) in the duct: 


f = (f ) 
c c mn 


— k 

2ir c 


- ij ¥ 2 + <?> 2 


(1/ sec) 


(50c) 


X = ' = ¥ * 
c 


2ab 


\J ( ? )2 + ( f )2 ^ (rob) 2 + (na)‘ 


(m) (50d) 


where c = f A 
c c 

There are a doubly infinite possible number of discrete acoustic wave 

modes corresponding to all the combinations of the integers m and n. 

From Eqn. (49) one sees that m represents the number of half-sine 

variations in the x-direction and n represents the number of the half-sine 

variations of the acoustic wave field components in the y-direction. 

Thus, the acoustic wave components of the mode (m,n) in the rectangular 

duct with rigid walls will be denoted by p and u and is given in 

mn mn 

Eqns . (49) . 

Of particular interest is the mode (m = 0, n = 0) which is referred 

. ■ } ■ . ■ : ■ 

td as the fundamental mode or the dominant mode of the acoustic wave in 

the rectangular duct with rigid walls. Substituting in Eqns. (50), 

m = 0 and n = 0, one obtains for the fundamental dominant mode: 


^C-'OO 0; 


(Voo - °; 


( Voo = " 


(51a) 
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From Gqns. (40), (41), (43) and (44), one obtains for m = 0 and n ■ 0: 


i (1) 

q o<T k = 


(V p>00 = ( VoO ~ c; 


( VoO " A; (z dl)0 ~ pC " Z o ; 


I. 

2 


(51b) 


Taking m = 0 n = 0 and qgQ= k in Eqn. (49), one obtains: 


p (x,y,z,t) = Re/2 p e i(kz-wt) - /2 p cos(kz-ait) 
00 0 ° 


(51c) 


U zoi X,y,Z,t) = Re ^ ^ P o e i(kZ " Ut) “ ^ ~ P 0 cos(kz-o)t) (51d) 




(51e) 


One finds that the fundamental dominant mode m = 0 and n = 0 in the duct 
with rigid walls given in Eqn. (51) is identical within the confined 
cross section of the duct with the acoustic uniform plane wave equations 
given for the infinite medium in Eqn. (23). Acoustic waves with every 
frequency can propagate in the duct in the fundamental mode (m = 0, 
n 0) since for all frequencies f > (f c )oo = This acoustic wave of 
the fundamental mode will propagate in the duct with infinite medium 
velocity, wavelength and impedance, as if the duct of the rigid walls 
were not there. The 0 = -j in Eqn. (51b) indicates that the fundamental 
dominant mode propagates parallel to the duct walls with no reflections. 

Let an acoustic uniform plane wave be propagating in an infinite 
fluid (air) in the z-direction. At each point in space it will have the 
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wave pressure scalar p, and the wave velocity vector u in the direction 

of propagation z. Both p and u will be constant along any phase plane 

z 

z = z q at any particular time t = t Q , and will not vary in x or y 

directions. Since u = 0 and u =0, one can introduce the four rigid 

x y 9 

walls of the duct of a rectangular cross-section directed in the z- 
direction without affecting the acoustic plane wave at all, since it 
does already obey the boundary conditions introduced by the rigid duct 
walls. Thus, the fundamental dominant mode of propagation of the acoustic 
wave in the duct with uniform rectangular cross section of rigid walls 
is that part of the same acoustic uniform plane wave in infinite medium 
which is confined by the duct walls without any changes in the pressure 
field, the velocity field, or any other parameters of propagation of the 
acoustic plane wave in an infinite fluid (air). 

Assuming the particular case where the cross section of the duct is 
a square a = b, then one has the cut off wave number for mode (m,n) 
from Eqn. (50a) in the form: 


( ^ c c ) m,n = \/ ( T )2 + ( f )2 = l\P + n2 (1/m) a=b (52a) 


where one has (k,) 0 Q = 0, (k^ Q = (k c > 0 1 = ~ and Eqn. (52a) may 
be rewritten in the form: 


^c^i»n / 2 

( Vi,0 = v ra 


+ n 


for a=b 


(52b) 


For the case of a square cross section a = b, the cut off frequency for mode 
(m,n) from Eqn. (50c) is in the form: 
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+ < f )J ' f.\P + " 2 *'•«> a - b < 52 '> 


where one has (f' c > 0 Q = 0, (f ( ,) 1 0 * ( f c > 0 l = 2a and Eqn ‘ ^ 52c ^ may 
be rewritten in the form: 


< f c> 


(f c > 


m,n [ 

i,0 v 


2 ^ 2 
m + n 


for a=b 


(52d) 


The numerical values of Eqn. (52b) and Eqn. (52d) are given in the 
Table on the following pages . 

If the signal frequency of the acoustic wave in the square duct 
c 

a = b is f < (-^ ~) , then only the fundamental dominant mode (0,0) 

can propagate and all the other modes will attenuate. If the signal 

frequency of the acoustic wave in the square duct a - b is (-^) < f < 1.4142(^0 

then only the acoustic modes (0,0) , (1,0) , (0,1) can propagate and all the 

other acoustic modes will attenuate. If the signal frequency of the 

c c 

acoustic wave in the square duct a = b is ^2a^ < ^ < 2.2361(— )» then only 

the acoustic modes (0,0), (1,0), (0,1) , (1,1) , (2,0) , (0,2) will propagate and 

all the other acoustic modes will attenuate. If the signal frequency of 

the acoustic wave in the souare duct a = b is f = 6.5(-r— ) then the 

Za 

acoustic modes (5, 4), (4, 5) and all the modes listed above them will 
propagate, since f > (f^)^ ^ for these modes, and the modes (6,3) (3,6) 
and all the modes listed below them in Table A will attenuate since 


f < (f ) for these modes, 
c m , n 
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TABLE A 


CUT OFF FREQUENCIES FOR MODES (m,n) IN SQUARE DUCTS 


( 0 , 0 ) 

( 1 , 0 ) 

( 1 , 1 ) 

( 2 , 0 ) 

( 2 , 1 ) 


mode (m,n) 


a = b 


(k ) 
c m,n 


oo 


c 1 , 0 


'c^m,n / 2 , 2 

v/m + n 

Vl,0 V 


( 0 , 1 ) 

( 0 , 2 ) 

( 1 , 2 ) 


REPRODUCIBILITY OF Tin 
ORIGINAL PAGE IS POOR 


0 

1.0000 

1.4142 

2.0000 

2.2361 


\ ' 

(2,2) 



j 2.8284 

■ 

(3,0) 

(0,3) 


! 3.0000 


(3,1) 

(1,3) 


• 

3.1623 

. . 

(3,2) 

(2,3) 


3.6056 

— 

(4,0) 

(0,4) 


4.0000 


(4,1) 

(1,4) 

1 


4,1231 


(3,3) 



4.2426 


(4,2) 

j 

(2,4) 


4.4721 


1(4,3) 

I 

(3,4) 

(5,0) 

(0,5) 5.0000 

1 

i i ■ ■ 

(5,1) 

(1,5) 


1 

j 5.0990 


;(5,2) 

(2,5) 


5.3852 

.(4,4) 



j 5.6569 

5(5,3) 

i : 

(3,5) 


■ 1 ■ ■ 

| 5.8310 

j(6,0) 

(0,6) 


1 

j 6.0000 

j( 6 , 1) 

(1,6) 


1 6.0828 

» - ■ ----- - 
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TABLE A ( cont inued ) 



mode 

(m,n) 

— 


_ i f c^m,n / 2 2 


HI 

( Vl,0 

.<V.i.o v 

(6,2) 

(2,6) 




6.3246 

(5,4) 

(4,5) 




6.4031 

'(6,3) 

(3,6) 




6.7082 

,(7,0) 

(0,7) 




7.0000 

,(5,5) 

(7,1) 

(1,7) 



7 . 0711 

(6,4) 

(4,6) 




7 . 2111 

(7,2) 

(2,7) . 


. 


7.2801 

(7,3) 

(3,7) 


• 


7.6158 

(6,5) 

(5,6) 




7.8102 

(8,0) 

(0,8) 




8.0000 

(7,4) 

(4,7) 

(8,1) 

(1,8) 


8.0623 

(8,2) 

(2,8) 




8.2462 

(6,6) 





8.4853 

(8,3) 

(3,8) 




8.5440 

(7,5) 

(5,7) 




8.6023 

(8,4) 

(4,8) 




8.9443 

(9,0) 

(0,9) 




9.0000 

(9,1) 

(1,9) 



: 

9.0554 

(7,6) 

(6,7) 

(9,2) 

(2,9) 


9.2195 

(8,5) 

(5,8) 




9.4340 









CHAPTER V 


VIBRATIONS OF THE PAN E L - BASIC THEORY 


In this chapter the equations governing the vibrations of a rec- 
tangular panel and its boundary conditions will be summarized as based 
on the theory of static and dynamic plates and the theory of elasticity. 

The detailed derivation of these equations may be found in two of the 
textbooks on the theory of plates listed in the Bibliography. In the 
first part the static plate will be considered and it will be extended later 
to dynamic analysis. 

The shape of a plate is defined by describing the geometry of its 
middle surface, which bisects the plate thickness "h" at each point. 

The small deflection plate theory developed by Kirchhoff and Love is 
based on the following assumptions: 

1. The material of the plate is elastic, homogeneous and isotropic. 

2. The plate is initially flat (in the x-y plane). 

3. h << a and h << b where "h" is the thickness of the plate (in the 
z-direction) and "a" and "b M are the dimensions of the plate 

(in the x-y plane) . 

4. p << h where rj is the lateral deflection or displacement 
of the plate (in the z-direction) and h is the thickness of 
the plate. 

5 . I-— -I << 1 ,|^| << 1 where ~ is the slone of the deflected 

1 Bx 1 1 By 1 Bx 


middle surface in the x-direction and is the slope of the 

By ■ 

deflected middle surface in the y-direction. 

The deformations due to transverse shear forces will be neglected. 
Straight lines, initially normal to the middle surface (in the 
z-direction) remain the same after the deformation. 
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6 \ 


7 




. The deflection of the plate is produced by displacement n of 
points of the middle surface normal to its initial plane 
(in the z-direction) . 

8. The stresses normal to the middle surface are neglected. 

9. The strains in the middle surface by inplane forces are 
neglected in comparison with the strains due to bending. 

The validity of the Kirchhoff-Love plate theory is assumed for the small 
deflections of the panel due to the acoustic waves incident on it. 

Taking the rectangular plate in the x-y plane, with the z axis 
perpendicular to it, the summation of all the forces in the z-direc- 
tion on a small element of area of the plate is zero, and it yields 
the first equilibrium equation: 


aq 3q 

x y 

8x 9y - P z 


(53a) 


where: p = the external load per unit area (or load density) (N/m ) dis- 

z 

tributed on the upper surface of the plate. 

= the transverse lateral shearing force per unit of length 
of the y-direction (N/m). 

- the transverse lateral shearing force per unit of length 
of the x-direction (N/m). 

Since the sum of all the moments acting on a small element of area 
around the y axis is zero, one obtains the second equilibrium 
equation: 


9m 9m 

+ — 2 * s 
9x 9y 


(53b) 
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where: m x = the bending moment around the y axis per unit length of 

the y-direction (N) 

m = the torsional twisting moment around the y axis per unit 
yx 

length of the x-direction (N) 

Since the sum of all the moments acting on a small element of area 
around the x axis is zero, one obtains the third equilibrium equation: 


dm dm 
— ! I + -SL = q 
dy ax y 


(53 c) 


where: m^ = the bending moment around the x axis per unit length of 

the x-direction (N) 

m = the torsional twisting moment around the x axis per unit 

xy 

length of the y-direction (N) 

Since there are no forces in the x and y directions and no moments 
with respect to the z axis, the three Equations (53) completely define 
the equilibrium of the element. Substituting Eqn. (53b) and Eqn. (53c) 
in Eqn. (53a) and taking m = m , one obtains: 


2 2 „ 2 
dm dm dm 

JL + 2 2Z. + J. 

dx 2 9x9 ^ dy 2 


By using Hooke’s law one may find: 


-D [^%+ v 

dx ay 


(55a) 


2 2 
nr 3n , lii 

m v = “ D t~2 + v 2 ] 

dy dx 


(55b) 


m = -(1 - v) D - - - 
yx dxdy 


(55c) 


where : D 


= Bending or flexural rigidity of the plate. (Nm). 

12(1 - O 2 

E = Young's modulus of elasticity (N/m ) which appears in Hooke's 
law. 

v = Poisson's ratio (v = 0.3 for steel and aluminum), 
h = thickness of the plate (m) . 

n = lateral deflection or displacement of the plate in the z- 
direction (m) . 

Substituting Eqn. (55) in Eqn. (54), one obtains: 


i_IL + 2 
8x 4 


.4 .4 p 

5 n , 3 n _ z 

2 2 + 4 D 

3x 3y 3y 


(56a) 


which may be also rewritten in the form: 


4 2 2 

dv n .« dv v n 


n , 9 . 3 ^ .. 

= D( — s- + — s-rri - 


(56b) 


3x 


8y 


Equation (56) was first obtained by Lagrange in 1811 and is called the 

Lagrange's equation of the non-homogeneous biharmonic equation. The 

static plate problem has been reduced to the solution of the linear 

partial differential equation (56) for the lateral deflection or 

displacement n = n(x,y), subject to boundary conditions of the plate. 

The transverse shearing forces q and q may be expressed in terms 

x y 

of the lateral deflections q(x,y) by substituting Eqn. (55) in Eqn. (53b) 
and Eqn. (53c): 


^x 


q y = 


2 2 

_ D r 9_n + UL] 

3 V [ 3 ,: 2 3 y 2 ’ 

_ D x riji + 

^ ‘ 3x 2 Sy 2 


. 3 „2 

■ D t; 7 n 


„ 3 „2 

-D t— V q 

3y 


(57a) 


(57b) 
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Substituting Eqn. (57) in Eqn. (53) , one obtains the differential equation 
(56). 

The Lagrange’s equation (56) is a fourth order partial differen- 
tial equation for the lateral deflection or displacement r\ = n(x,y). 

In order to obtain a unique solution, one should give two boundary 
conditions for each boundary or edge of the plate at x = 0, x = a, y - 0 
and y = b, depending on the kind of support the edge has. Some different 
possible boundary conditions at the edge for the different cases are 
discussed in the following: 

A. Simply supported edge - If the edge x = a of the plate is 

pimply supported, the deflection n along this edge must be zero. At 

I 

the same time, this edge can rotate freely with respect to the edge 
line, i.e., there are no bending moments m^ along this edge. The 
boundary conditions for this edge are, using Eqn. (55a): 


n = 0 and m 
u x 


-D[^.+ v *-%] = 0 
3x 3y 


at x = a (58a) 


Since at x = a one has n = 0, and therefore -r^- = 0 and — ^ = 0 vanish 

d ' 3y 

as well, one may rewrite Eqn. (58a) alternatively: 


n - 0 and =0 at x = a (58b) 

ax 

Similarly, one has from Eqn. (55b) for the edge y = b: 

2 2 

q = 0 and m = -D[-Mr + v — x] =0 at y = b (58c) 

7 3y^ 3x 


or alternatively in the form: 


n = 0 and ^—9- = 0 at y = b . (58d) 

ay 

Equations (58) represent the boundary conditions for a simply supported 
edge. 

B. Clamped edge (also called fixed edge or built-in edge) - In 
this case the deflection along this edge is zero, and the tangent 
plane to the deflected plate along this edge coincides with the initial 
position of the plate. Therefore the boundary conditions for the 
clamped edges are given by: 

tv = 0 , = 0 at x = a (59a) 

oX 

P “ 0 , |y = 0 at y = b (59b) 

C. Free edge (also called the statical boundary condition) - In 
this case the edge moment and the transverse shear force (v'l are 
zero. The shearing force at the edge of the plate consists of the 
transverse shear and the effect of the torsional twisting moment. The 
boundary conditions at a free edge may be given in the form: 

8m 

m =0 and v = q + =0 at x = a (60a) 

x x x 8y 

Substituting Eqns. (55a), (55c) and (57a), one obtains the boundary 
condition for the free edge in terms of lateral displacement n as follows 

2 2 3 3 

m = -D [— — r + v — x] = 0 and v = -D[-~ ?• + (2-v) — - j ] = 0 
x 8x 8y x 3x 8x9y 

at x = a (60b) 
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Similarly, one has for the edge at y = b: 


3m 

m = 0 and v = q + — = 0 

y y y 9x 


m 


' . 2 * 2 

-D[^ + v ~ S'] = 0 and v 


-b[*4 


3y 


ax 


3y- 


at y - b (60c) 


+ (2 - v)-^-^]= 0 
SyBx 

at y = b (60d) 


Equations (60) represent the boundary conditions for a free edge. 

In the dynamic analysis of the vibrations of the plate only the 
literal motion in the z-direction is of interest, since the rotational 
inertia effects may be neglected. Therefore only the inertial forces, 
associated with the lateral translation of the plate in the z-direction 
will be considered. The partial differential equation governing the 
vibrations of the plate may be derived by applying the d’Alembert’s 
dynamic equilibrium principle. The inertia force density associated 
with the lateral translation Tl(x,y) of a plate element can be expressed 
by: 


2 

inertial , d n 

p = -p h — ^ 

* P 3t 2 


(61) 


where: p ■*- nertia ^ = the inertial force per unit area (N/m^) 

z 

n = lateral displacement of the plate (m) 

3 

p = mass density of the plate material (Kg/m ) 

P 

h - thickness of the plate (m) 

2 

pph - mass density of the plate per unit area (Kg/m ) 
Extending the differential equation of static equilibrium (56b; 
by adding the inertial force Eqn. (61) , one has : 
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tv* _ „ . n inertial . 3 rt 

°v n - p z + p z “ p z - Pp h — i 

F 3t 


which can be rewritten in the form: 


<x.y,t> + p h . p ( X>y ,t.: 

P 3t Z 2 


tO (62) 


where Eqn. (62) represents the inhomogeneous partial differential - 

equation ofn forced, undamped motion of the plate. 

For the case of a freely vibrating plate, the external applied 

force p (x,y,t) = 0, and Eqn. (62) becomes: 
z 

DV 4 n(x,y,t) + , h . 0 (63) 

p at 

where Eqn. (63) represents the homogeneous partial differential equa- 
tion of free, undamped flexural vibrations of the plate. Assuming 
harmonic time variation of the lateral displacement and of the external 
force per unit plate area: 


n(x,y,t) = 


Re/2 n(x,y) 


-ioj't 

e 


(64a) 


p (x,y,t) = Re/2 p (x,y)e" 1U3t (64b) 

Z ^ Cd 

and substituting in Eqn. (62) and Eqn. (63) , one obtains: 

4 2 

DV n(x,y) - p hoi n(x,y) = p(x,y) (65a) 

~ P ~ fZ 

4 2 

DV n(x,y) - p htD r,(x,v) = 0 


(65b) 



where q(x,y) and p (x,y) are complex phasors of the lateral displace- 
^ z 

ment and the external force per unit plate area; the complex sign 
will be omitted when it is obvious . 

The natural resonance frequencies of the free, undamped flexural 
vibrations of a simply supported rectangular plate of dimensions "a" and 
"b" may be found by assuming for the lateral displacement q of the 
homogeneous equation (63) the solution originally proposed by Navier 
in 1820 for the static plate: 


n(x,y,t) = Csin sin ^’ a>t (66a) 

where m and n are integers. Eq. (64a) obeys the boundary conditions 
for a simply supported rectangular plate at all four edges, as given 
in Eqn. (58b) and Eqn. (58d) . Substituting Eqn. (66a) in Eqn. (63), one 
obtains: 


(=> 4 + 2 (= J 2 (- f ) 2 + (^) 4 - u , 2 = 

= [(— ) 2 + C ^) 2 ] 2 - ^ 0) 2 = o ( 66 b) 

a b D 


where 03 = 00 = 2irf is the resonance frequency of the mode (m,n) of 

mn mn 

the plate vibrations. From Eqn. (66b) , one obtains the resonance fre- 
quencies of the simply supported rectangular plate as follows: 


00 

mn 


. mTT v 2 

= 


(f) 2 ] 



f 

mn 


= f[(f ) 2 + 


(£) 2 ] 



(1/sec) (66c) 

(1/sec) (66d) 
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where: a,b ■ dimensions of the plate (m) 

h * thickness of the plate (m) 

3 

pp = mass density of the plate material (Kg/m ) 

D = bending or flexural rigidity of the plate (Nm) 

The shape of the different modes of vibrations of the plate are given by 
Eqn. (66a) for the simply supported plate, while the corresponding 
resonance frequencies are given by Eqn. (66c) or Eqn. (66d) . The fun- 
damental mode of flexural vibration has a single sine wave in the x-direc- 

tion and a single sine wave in the y-direction. The pertinent resonance 
£ 

frequency f- .. may be found by taking m = 1 and n = 1 in Eqn. (66d) . 

1 » X 

For a simply supported square plate where a = b the mode (m = 2, 
n » 1) and the mode (m = 1, n = 2) have different shapes, but the same 
resonance frequency f.. 9 * f - . ; such modes of flexural vibrations are 
called degenerate inodes . For a rectangular plate a ? b one has degenerate 
modes if the (a/b) ratio- is a rational number. 

Let us consider a simply supported rectangular plate of dimen- 
sions "a" and "b," acted upon by a given external arbitrary lateral load 

density p (x,y), which varies harmonically in time with circular fre- 
z 

quency u>... One has, therefore, using complex phasors: 

P z Cx,y,t) = p z (x,y)e~ 1U3t (67a) 

where p (x,y) is a given function and can be expanded in a double Fourier 
z 

series of sine terms as follows : 

P (x,y) - Z Z P sin sin (67b) 

z , , m,n a b 
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(67c) 


a b 

p = — - f [ p (x,y) sin sin ^^d xdy 
m,n ab f / ■ z' * J/ a b 

0 0 

In accordance with the method of Navier, the particular solution for the 
Lagrange equation in Eqn. (62) for the simply supported plate is given 
in the form, for the harmonic time varying case: 


n (x,y,t) 


Z 

m=l 


Z 

n=l 


A sin 
m,n 


mrrx 

a 


sin 


SIZ 

b 


-iwt 

e 


( 68 ) 


where A are the unknown constants to be found for the particular 
m,n 

solution. Substituting Eqn. (67) and Eqn. (68) in the Lagrange equation 
(62) and matching coefficients, one obtains: 


A. _ [ D ( 


2 2 
m ir 


2 2 
n it 


m,n 


) 2 - p hoi 2 ] = P 

p ra,n 


(69a) 


where A may be found from Eqn. (69a) while P is calculated from 
m,n J m,n 

Eqn. (67c). Substituting Eqn. (66c) in Eqn. (69a), one obtains: 


p h A [*d 

K p m,n m,n 



(69b) 


or in the form: 



P 

m,n 


p h(w 
K p m,n 



(69c) 


where A in Eqn. (69c) is expressed in terms of the resonant frequency 
m,n 

u) and the load frequency u), to be substituted in Eqn. (68) for 
m,n 

the final particular solution in the form of a Fourier series double sum. 
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I 

l 

The homogeneous solution of the homogeneous equation given in Eqh. 

(66a) should be added to the particular solution of the inhomogeneous 
equation given in Eqn. (68) and Eqn.(69), to give the general solution 
subject to given initial time conditions* 





CHAPTER VI 


ACOUSTIC PLANE WAVE INCIDENT ON AN INFINITE PANEL 


In the present chapter the case of a plane acoustic wave in an 
infinite fluid (air) , normally incident on an infinite panel (plate) 
will be discussed. 

Assuming harmonic time variation e la)t , the inhomogeneous partial 
differential equation which governs the lateral displacement of the 


panel (plate) given in Eqn. (65a) can be rewritten as follows: 


v n - 


(70a) 


where: q(x,y) = lateral displacement of the plate. (m) 

p (x,y) = external net force per unit area in the positive 
z 

2 

z-direction. (N/m ) 

3 

pp = mass density of the plate. (Kg/m ) 
h = thickness of the plate. (m) 

gj = 2nf = circular frequency of the wave. (l/sec) 

Eh 3 

D = — = bending or flexural rigidity of the plate. 


12(l-v ) 


E = Young’s modulus of elasticity. (N/m ) 
v = Poisson’s ratio (v = 0.3 for steel and aluminum) 


Equation (70a) may be rewritten in the form: 


v n - y n = 5 p z 


(70b) 


where the plate wave number y is defined by: 


. p ha 2 12(l-v 2 )p co 2 

Y 4 = _E _E — 

" ° Eh 2 


(1/m^) 


(70c) 


Assuming in Eqn. (70b) no external force density p “0 and variation in 
only the y-direction n = n(y), one will obtain: 

^4 - y\ - 0 (71a) 

dy 

which will be solved to give the following free motion elastic plate 
waves of the lateral displacement of the plate propagation along the 
plate: 

-3 

n = n Q e ±Yy la)t attenuation plate elastic waves (71b) 

n = q propagating plate elastic waves (71c) 

where for tne propagating wave Eqn. (71c) the plate wave number can be 
given in the form: 

2 

- f = ¥ ■ t - 22(1 : >Pp ^ (i/ "> <72a > 

P P P 

where: c = wave velocity of the elastic plate waves. (m/sec) 

P 

X = wavelength of the elastic plate waves. (m) 

P 

£ '= frequency of the wave. (1/sec) 

From Eqn. (72a), one obtains: 

c - [ \ fiw = [ r-rs ]* /2irhf (m/sec) (72b) 

P 12(l-v )p 12(l-v ) p 

P P 
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A = [- 


12(l-v )p. 


i % ■ 


[“ 


12(l-v )p, 


On) 


(72c) 


where c = i\ and c = c (f). The plate is dispersive medium for the 
P P P P ' 

transverse displacement waves propagating along the plate. 

Let an infinite panel (plate) be situated at z = 0 in the x-y plane. 
Let an acoustic plane wave propagating in the positive z-direction be 
normally incident on the infinite plate and be given by Eqns. (23) in the 
form: 


i(kz-u)t). 

r I e 


U ; 
Z1 



z < 0 


(73a) 


This incident acoustic plane wave will be reflected by the infinite 
plate in the form of a reflected acoustic plane wave propagating in the 
negative z-direction and given by Eqns. (24) in the form: 


P = 

*r R 


-i(kz-f-o)t) 


u 


zr 


r 

pc 


z < 0 


(73b) 


The incident acoustic plane wave will cause the infinite plate to 
vibrate harmonically in the lateral positive z-direction in the form: 



z = 0 


(74a) 


where n is the lateral displacement of the plate in the positive 
z-direction ; the velocity of the plate in the positive z-direction 
may be found from Eqn. (74a) to give: 



dn 

dt 


-iwAe 


-icot 



(74b) 


No plate boundary conditions are applied to Eqn. (74a) since the plate 
is of infinite dimensions in both x- and y-directions. 

The vibrations of the plate given in Eqn. (74) will generate on 
the other side of the plate, z > 0, an acoustic transmitted plane wave 
in the positive z-direction which will be given by Eqns. (23) in the 
form: 



i(kz-a)t). 

x e 


u 


z t 



z > 0 (75) 


The total external force per unit area p z on the plate at z = 0 in 
the positive z-direction by the incident acoustic wave Eqn. (73a), the 
reflected acoustic wave Eqn. (73b) and the transmitted acoustic 
wave Eqn. (75), is given by: 


_ -ioit -liot 

Pz = P i e +P R e 


- 

T 


-iwt 


(76a) 


where the incident and reflected acoustic waves at z < 0 pressure the 

plate in the positive z-direction and the transmitted acoustic wave at 

z > 0 pressures the plate in the negative z-direction. Substituting 

4 

Eqn. (74a) and Eqn. (76a) in Eqn, (70b), one finds V n = 0 and it 
becomes : 


-A = i ( P I + P R - P T> 


(76b) 


where the factor e 1Wt has been cancelled, 
may rewrite Eqn. (76b) in the form: 


4 2 

Since y D = P^hoj , 


one 
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(76c) 


'r-zze.x ....... ... 




V“ a-p i + p r -p t 


The plate velocity u and the fluid (air) velocity of the acoustic 

zp 

waves on either side of the plate should be identical. On the positive 

side of the plate z > 0 at the plate z = 0, denoted by z = 0 + , the 

plate velocity u should be identical with the transmitted acoustic 
zp 

wave velocity vector n as follows: 


u z p (z = V = u zt (z = V 


(77a) 


Substituting Eqn. (74b) and Eqn. (75) in Eqn. (77a), one obtains: 


-iajA = P 
pc I 


(77b) 


On the negative side of the plate z < 0 at the plate z - 0, denoted by 
z = 0 , the plate velocity u should be identical with the sum of the 
incident acoustic wave velocity vector u ^ and the reflected acoustic 
wave velocity vector u ;; as follows: 


u (z “ 0 ) - u . (z = 0 ) + u (z - 0 ) 
zp - zi - zr - 


(78a) 


Substituting Eqns. (73a) , (73b) and (74b) in Eqn, (78a), one obtains: 


■iwA = — P T - — P„ = — (P T - P B ) 
pc I pc R pc I R 


(78b) 


From Eqn. (77b) and Eqn. (78b), one obtains: 
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p = p + p 
I R T 


( 79 ) 


The pressure of the incident wave is divided among the pressure of 
the reflected wave and the pressure of the transmitted wave. 

From a given acoustic incident plane wave P^, on the infinite 
plate, one may obtain the unknown P^, P T and A by solving the three 
linear inhomogeneous equations (76c), (77b) and (78b). 

From Eqn. (77b), one obtains: 


P T = -impc A 


(80a) 


By substituting the additional equation (79) in Eqn. (76c), one 
obtains : 


p hu 

P = - — 2 A 

R 2 


(80b) 


Substituting Eqns. (80) in Eqn. (79), one obtains: 


p hoi 


p ha) 


Pj = C-iitfpe — > A = ~iwpc(l - i -^~)A 


(80c) 


The coupling parameter y between the plate and the air is defined 
as follows : 


P = 


2p 

p P h 


(1/m) 


where: p = air density. (Kg/m ) 

pp = plate material density, 
h = thickness of plate. (m) 


(Kg/m ) 
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By substituting the parameter y = 

2 Pp 
k = — « — - in Eqn. (80c), one obtains: 

C A 


and the acoustic wave number 


P T = -iajpc (1 - i — ) A 

y 


(80d) 


where (— ) is a non-dimensional number. From Eqn. (79) and Eqn. (80), one 
obtains: 


* _ iHjP'C I 

A " . k 

1 - 1 y 


P T 


P R = 


1 * k 
1 "Vy 

-V' P i 

1 - 1 y 


(81a) 


(81b) 

(81c) 


The Transmission Loss coefficient is defined as the ratio of the 
incident acoustic wave pressure power to the transmitted acoustic wave 
pressure power. From Eqns. (26), (73a) and (75), one obtains the 
Transmission Loss coefficient in decibels in the form: 


TL db - 10 log 


= 10 log 


(82a) 


where TL,, is the Transmission Loss coefficient expressed in decibels, 
db 

Substituting Eqn. (81b) for the present case in Eqn. (82a), one has: 


TL db - 10 log |1 - 1 ^| 2 


10 log [1 + (£) 2 ] 


= 10 log [1 + ( 


up h , 


(82b) 
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Usually the measurements are done by two microphones which measure 
the pressure power of the acoustic waves. One microphone is situated 
on the source side of the panel at z = -d^ (the source microphone) 
and it measures the total acoustic wave pressure power at that point, 
bdth the incident and the reflected acoustic waves. The other micro- 
phone is situated on the other side of the panel at z = +&2 (the 
receiver microphone) and it measures the total acoustic wave pressure 
power at that point of the transmitted acoustic wave. The Noise Reduc- 
tion coefficient is defined as the ratio of the total acoustic wave 
pressure power measured by the source microphone at z = -d^ to the acous- 
tic wave pressure power, measured by the receiver microphone at z = 

From Eqns . (73) and (75), one obtains the Noise Reduction coefficient 
in decibels in the form: 


NR^ = 10 log 


P i l ’ P r z=-d. 


l p tl 2 =+d. 


= 10 log 


P T e‘ ikd l + P,,e +ikdl | 2 

1 K 


V 


ikd2 


(83a) 


where NR., is the Noise Reduction coefficient expressed in decibels, 
db 

Substituting Eqn. (81b) and Eqn. (81c) for the present case in Eqn. (83a) 
one obtains: 

i k 


NR,, - 10 log 
db 


-ikd'i u ^+ikdi 


i-i * 

L_ 


1-i * 
y 


-•=10 log 


.k, -ikdi .k +ikdi 
(1-x— )e A -i— e i 

y y 


= 10 log |e - i — 2 cos kd | 

y 1 


. 2k 


=10 log coskd, -isin kd., - i — cos kd- 
1 J- u - 1 - 
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2 2k 2 

10 log [cos kd- + (sin kd- + — cos k dj ] = 
1 1 u 1 


2k 2 2 

10 log [1 + — sin kd- cos kd.. + ( — ) cos kd ] = v 
u 1 l u i 


10 log [1 + — sin 2 kd n + (— ) 2 cos 2 kdj 
u 1 u 1 


(83b) 


where the Noise Reduction coefficient NR^ b depends on the position of 

the source microphone z = -d . For the particular case kd.. = 
d l 

2u — << 1, Eqn. (83b) will reduce to: 


?v ■? up i, h 7 

NR db = 10 log [1 + (~ ) Z ] = 10 log[l + (-^~) 1 


for |kd 1 |« 1 


(83c) 


for the case of kd. << 1 the Noise Reduction coefficient NR,, does 

1 db 

not depend on the position of the source microphone z = 

The Insertion Loss coefficient is defined as the ratio of the 
incident acoustic wave pressure power measured by the receiver micro 
phone without the panel, to the transmitted acoustic wave pressure 
power measured at the same point by the receiver microphone with the 
panel installed. From Eqns. (73a) and (75), one obtains the Inser- 
tion Loss coefficient in decibels for the present case in the form: 


!L db = 10 log 


= 10 log 


(84a) 


where IL,, is the Insertion Loss coefficient expressed in decibels, 
db 
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Substituting Eqn. (81b) for the present case in Eqn. (84a), one has: 


If 2 lr 9 o 

IL db = 10 log |1 - i = 10 log[l + (^H = 10 log[l +(~ 2 ^-) ] (84b) 
From Eqns. (82b) and (84b), one finds for the present case: 

"■db ■ IL db (85a) 

From Eqns,. (82b) and (83c), one finds that for (k/y) >> 1 one has: 

( 2k) 2 

NR . - TL = 10 log A 2 - 10 log(£) 2 = 10 log = 

uu ^ ^ V o 

V 

= 10 log(2 2 ) = 20 log 2 = 20 X 0.3010 = 6 db 

which may be rewritten in the form: 

NR db = TL db + 6db = IL db + 6db for k>> y and | kdj «1 (85b) 

Equation (82b) is generally known as the "normal incidence mass law," 

For the case k >> y Eqn, (82b) becomes: 

, oop h TTp hf 

TL db = 20 log (-) = 20 log = 20 log (— |— ! ) for k»y (86a) 

From Eqn, (86a), one sees that when the plate surface density (p^h) is 
doubled, the Transmission Loss coefficient increases by 6 db. One 
also sees that if the frequency f doubles (= the frequency is raised by 



one octave) the Transmission Loss coefficient increases by 6 dh. 

Eqn. (86a) may be rewritten in the form; 

it p h 

TL ,, = 20 log ( — £-) + 20 log f = C + 20 log t for k» y (86b) 
db pc 

Thus if the Transmission Loss coefficient in decibels is drawn against 
(log f ) , it will give a straight line, with a slope of 20 (= doubling 

the frequency will increase TL^ by 6 db) and it will have the value 
TL^k = C at the frequency f = 1 (1/sec), or the value TL^ = (C+40) 

at the frequency f = 100 (1/sec) . In accordance with Eqn. (85b) the 

Noise Reduction coefficients NR,, for k » p will be a straight line 

db 

parallel to the TL^ straight line, and 6 db above it, For the case 
k << p one finds from Eqns. (82b) and (83c) that both the TL^ 

coefficient and the NR,, coefficient should approach the value zero, 

db 

and the normally incident acoustic plane wave for very small frequencies 
will propagate through the infinite panel as if it was completely 
transparent to it. 
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CHAPTER VII 


ACOUSTIC PLANE WAVE INCIDENT ON A FINITE PANEL 

In the present chapter the case of a plane acoustic wave normally 
incident on a finite clamped panel (plate) in a rigid duct will be 
discussed. 

It has been shown before that the fundamental mode of propagation 
of the acoustic waves in a rigid duct are identical with the plane 
acoustic wave propagating in the infinite fluid (air) limited by the walls 
of the rigid duct. The introduction of the walls of the rigid duct 
parallel to the direction of propagation does not affect the plane 
acoustic wave, since the boundary conditions of the walls of the rigid 
duct are obeyed; the fundamental mode of propagation of the acoustic 
wave in the duct propagates in the duct as a plane acoustic wave, as if 
the duct rigid walls were not there. Thus, the introduction of the 

duct rigid walls parallel to the direction of propagation does not 

alter the acoustic waves discussed in the previous chapter. 

When one introduces a finite rectangular clamped panel (plate) in 
the duct, one should apply the boundary conditions for a clamped 
rectangular plate as discussed previously. One is not able to use the 
solution for the infinite panel (plate) discussed in the previous 
chapter without any boundary conditions. Thus, the introduction of the 
duct rigid walls parallel to the direction of propagation does not 

affect the plane acoustic waves discussed in the previous chapter, but 

it does introduce the effect of the clamped boundary conditions of the 
panel (plate), which is now clamped to the rigid wall of the duct. The 
solution given in the previous chapter for the infinite panel (plate) 


59 


for the lateral displacement n in Eqns. (74a) and (81a) do not obey 
the boundary conditions for the clamped edges of the panel (plate) at 
x = 0, x = a, y = 0, y * b as given in Eqn. (59). Therefore the 
solution for the infinite panel (plate) given in the previous chapter 
should be altered. 

In order to simplify the analysis and limit it only to the funda- 
mental resonance frequency of the panel (plate) in the present chapter, 
it is assumed that the effect of the four clamped edges of the plate is 
equivalent to a spring, forcing the plate to return to its position of 
equilibrium without acoustic plane wave incident on it. Thus, the 
plate has been idealized by introducing the effect of the clamped edge 
boundary condition as an equivalent system of a single degree of freedom 
only. For small lateral displacements of the plate n, the spring-like 
force per unit area is proportional to the lateral displacement of the 
plate and can be expressed in the following form by using the stiffness 
constant of the plate: 


spring 

z 



(87) 


where: p s P r *- n § = spring like force density. (N/m^) 

z 

q = the lateral displacement of the plate. (m) 

3 

K = the stiffness constant of the plate. (N/m ) 


The formulation of the structural damping of the vibrating plate can 
be accomplished by expressing the damping force per unit area as a 
purely imaginary constant , being proportional to the small lateral 
displacement, as follows, where i = /-l. 
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iotKn 


( 88 ) 


damping 


where: p dain P in S = the damping force density. (N/m 2 ) 

z 

n = the lateral displacement of the plate. (m) 

3 

K = the stiffness constant of the plate. (N/m ) 
a = the damping factor. (non dimensional) 

Equation (88) has been employed extensively in the dynamic analysis 
of aerospace structures. 

The inertia force density associated with the lateral translation 
q(x,y) of a plate element has been given in Eqn. (61) as follows: 


P 


inertial 

z 



(89) 


where: p inertia] - = the inertial force per unit area. (N/m 2 ) 

z 

3 

Pp = mass density of the plate material. (Kg/m ) 
h = thickness of the plate. (m) 

2 

Pph = mass of plate per unit area. CKg/m ) 

Using Eqns. (87) and (89), the natural fundamental circular resonance 
frequency will be defined for a spring-like system of one degree 
of freedom of this type as follows: 



Extending the differential equation of static equilibrium, Eqn. (56b), by 
adding the spring force, Eqn. (87), the damping force, Eqn. (88), and 
the inertial force, Eqn. (89), one obtains: 
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-.4 , spring . damping , inertial 

DV n ■ P + P„ + P, + P, 


z z 


3 2 n 


- P z - K n - ictKn - p h ^ 


3t 


which can be rewritten in the form: 


DV^n + p^h + K(1 + ia)n - P 2 (x»y,t) 


(91a) 


3t 


Substituting from Eqn. (90) K = p^hu^ " in Eqn. (91a), one has: 


DV^n + P_h[-^— 5- + a) 2 (1 + ia)n] = P_(x,y,t) 
P 3t ° 


(91b) 


where Eqns.(91) represent the inhomogeneous partial differential 
equation of forced and damped motion of the clamped plate, taking only 
the fundamental circular resonance frequency a) Q into account. 

Assuming harmonic time variation e as in Eqn. (64), one 
will obtain from Eqn. (91b) : 


4 2 00 

DV n - p htu [1 %r- (1 + ia)]n = p (x,y) 

p z z 

00 


(92a) 


which may be rewritten in the form: 


4 4 1 

v n - Y o n = 5 p z 


(92b) 


where the plate wave number y^ for the present case is defined by: 


■2 2 2 
, p ha) oo / oo / 

Y = -V- [1 - -V (1 + ia)] - Y [1 - ~ (1 + ia) ] (1/m 4 ) (92c) 


D 


oo 
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p hu) 2 12(1 - V 2 )p W 2 


(1/m 4 ) 


oj = circular frequency of th& external force. (1/sec) 

in = fundamental circular resonance frequency of the plate. (l/sec) 
o 

a = the damping factor. (non-dimensional) 

Eh 3 

D = — = bending or flexural rigidity of the plate. (Nm) 

12 (1-v ) 

For the particular case K = 0, one has a) = 0 and = Y* 

o o 

Let the clamped panel (plate) be situated at z = 0 in the x-y 
plane in the rigid duct. Let the fundamental or dominant mode 
(m = 0, n = 0) of the acoustic wave propagating in the positive z- 
direction in the rigid duct be normally incident on the clamped plate and 
be given in the form: 


p i = p i e 


i(kz-iot) 


P i 

u = — 
zx pc 


z < 0 


(93a) 


The reflected fundamental mode from the clamped plate will propagate 
in the negative z-direction and is given by: 


P = p D e 
r R 


-i(kz-hi)t) 


z < 0 


(93b) 


The incident acoustic wave will cause the clamped plate to vibrate 
harmonically in the lateral positive z-direction in the form: 


n = Ae 


dn . . -io)t 
u = — = -ltuAe 
zp dt 


z = 0 


(93 c) 
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The vibrations of the plate in Eqn. (93c) will generate on the other 
side of the plate z > 0 the fundamental mode of the acoustic transmitted 
wave in the positive z-direction in the form: 


. i(kz- u t) = _t 

p t T * zt pc 


z > 0 (93d) 


Substituting Eqns. (93) in Eqn. (92b), one finds similarly to 


Eqn. (76b) : 


- >1 + P R - V 


(94a) 


Similarly to Eqn. (77b) and Eqn. (78b), one also finds from Eqns. (93) 


-io>A = — P T 
pc T 


(94b) 


■*“* ’ ^ (P I ' P R> 


(94c) 


From Eqns. (94b) and (94c), one also obtains: 


P = P + P 
I R T 


From Eqn. (94b), one obtains: 


P T = -ioipcA 


Substituting Eqn. (95) in Eqn. (94a) , one obtains: 
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DY 


P R “ " 


(96b) 


Substituting Eqns. (96) in Eqn. (95), one obtains: 


DY 


DY 


P I = (- iUpc ‘ * -^P^ 1 - 1 2 o^c )A 


(96c) 


From Eqns. (96), one obtains: 


A = 


— P T 

upc I 


1 - i 


°Y, 


2copc 


(97a) 


= 


1 - i 


Dy 


2copc 


(97b) 


Dy 


-x 


P R = 


2 cope I 


1 - i 


Dy, 


2copc 


(97c) 


The Transmission Loss coefficient defined in Eqn. (82a) may be 
found from Eqn. (97b) to give for the present case: 


TL db - 10 log 


= 10 log 


1 - i 


Dy 


4 2 


2copc 


(98a) 


Substituting Eqn. (92c) in Eqn. (98a) and defining k = — and y = 


one has: 
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U) 2 

TL .. = 10 log 1 - i -[1 - Ar (1 + ia) 1 
do ii * 


2 2 

t_ o) a t w *■ 

- io log la---^- ) - i - a - — )| ■ 

M aj a) 

2 2 
k 2 k 2 2 

= io log [a - ) + (-) a - -%-) i 

11 «, 2 “ » 2 


For the particular case of an infinite plate K = 0, therefore u> o 

and Eqn. (98b) becomes identical with Eqn. (82b) . 

The Noise Reduction coefficient defined in Eqn. (83a) may be 
found from Eqn. (97) to give for the present case: 


^ + Pr , 


^-ikdi +ikdi 


10 

log 

(1 

10 

log 

i -i 
e 


2 cos kd 






2 p 


Substituting Eqn. (92c) in Eqn. (99a) and defining k = — and p = — b 
Eqn. (99a) may be developed as follows: 


NR db = 10 log 


i ^k^l _ ^ 2k ^ _ __2_ ^ + ia)] cos kd- 
u J- 


0) 


10 log 


= 10 log 


?V u 0 W a 

cos kd- - i sin kd^ cos kd^[i(l 2 ~) + ] 


= 10 


2v w n a 2k o 

cos kd.. [1 r- ] - i[sin kd 1 + — (1 5 ~)cos kd.. ] 

1 ^ a ) 2 1 v a ) 2 1 

f 2 2k M n a 2 2k U o 2 ) 

log t cos [ 1 y 3 + [sin kd^ + — Hi 2 "^ cos ^d^] r(99b) 

u a) x y 0 ) J 


For the particular case of an infinite plate K = 0, therefore a> - 0, 
Eqn. (99b) becomes identical with Eqn. (83b). 

d l 

For the particular case that kd^ = 27r-^— << 1, one has cos kd^ = 1 
and sin kd^ <= 0 and Eqn. (99b) becomes: 


2 2 
/}•» Cli 0i ry ryi 2 (0 O 

NR db = 10 log [ (1 - — ) 2 + (— ) (1 - -y) 3 


where |kd^ << 1 


(99c) 


For the particular case of an infinite plate K = 0, therefore oj q = 0, 

Eqn. (99c) becomes identical with Eqn. (83c). For the case of 
|kd^| << 1 given in Eqn. (99c), the Noise Reduction coefficient NR db 
does not depend on the position of the source microphone z = -d^. 

The Insertion Loss coefficient in decibels defined in Eqn. (84a) 
may be found from Eqn. (97b) to give for the present case as in Eqn. (98b): 
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2 


r . 

p. 

2 

P T 

2 

DY 0 4 

11 IL db - 10 l °s 

p t 

= 10 log 

I 

P T 

= 10 log 

^ ^ 2o3pC 


k^a a 

10 log [ (1 - ~ "V 

03 


) 2 + (— ) 2 (i - 


03 


f) 2 i 


(100) 


where IL,, is the Insertion Loss coefficient expressed in decibels, 
do 

From Eqn. (98) and Eqn. (100), one finds for the present case: 


TL db " IL db 


From Eqn. (98b) and Eqn. (99c), one finds that for: 


(101a) 


i u> 

— (1 — ) 

U U 2 ' 

to 


> > 


k V 

y 2 

0) 


one has: 


NR db - TL db = 10 log 4 = 20 log 2 • 6 db 


which may be rewritten in the form: 


NR,, = TL ,, + 6 db = IL,, + 6 db 
db db db 


(101b) 


Assuming that in the region of frequencies oj considered in this 


report |kd | << 1 and substituting k = — in Eqn. (99c), one obtains: 
1 c 


2 2 ' 2 
O 0) 01 9 9 9 03 03 9 

NR,, = 10 log [ (1 - — — ) 2 + (— ) 2 ( ' 2 

db ° lie 03 yc oj 


) ] 


(102a) 


Since oj = 2irf and w = 2irf , one may rewrite Eqn. (102a) in the form: 

o o 


f z a 


f - _ f 2 


NR db = 10 log [(1 - £ _ )2 + (^) 2 (— ) 2 ] (102b) 
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where : NR^ 

u = 

f = 

u. = 
o 

f = 
o 

c = 



= Noise Reduction coefficient in decibels. 

the circular frequency of the acoustic wave. (1/sec) 

the frequency of the acoustic wave. (1/sec) 

the circular fundamental resonance frequency of the plate. 

(1/sec) 

the fundamental resonance frequency of the plate. (1/sec) 
the velocity of the plane acoustic wave in infinite 
medium. (m/sec) 

the damping factor of the vibrating plate. (non dimensional) 
the coupling coefficient y. (1/m) 


The Noise Reduction curve in Eqn. (102b) as a function of fre- 
quency f can be divided into three major regions: 

A. The mass controlled region : This is the region where the 

frequency of the acoustic wave f is much larger than the fundamental 

2 2 

resonance frequency of the plate f Q such that f >> f Q . In this 
region Eqn. (102b) becomes: 


m db ■ 10 lo * l<£> 2 f2 l - 20 ’•OS [ft £ 1 - 

= 20 log (— ) + 20 log f 

yc 

2itp h 

= 20 log ( 2-) + 20 log f 

PC 

= C + 20 log f for f 2 » f^ 2 (103) 
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In this region doubling themass per unit area of the plate (pph) will 
increase NR^ by 6 decibels, and raising the frequency f of the acoustic 
wave by one octave (= doubling the frequency f) will increase NR,, 

db 

by 6 decibels. Drawing NR^(log f) on a log-log paper will give a 

straight line with a slope of 20 (increase of NR,, by 6 db for every 

db 

octave of f) and a constant C such that: 

2irp h 

™db< f ■ » ’ C - 20 * - 20 108 <2^ 

p 

NR db (f = 100) = C + 40 db 

B. The damping controlled region : This is the region where the fre- 

quency of the acoustic wave f is approximately the same as the funda- 
mental resonance frequency of the plate f such that f = f . Since 
| a | << 1 and the first term in Eqn. (102b) varies slowly with f as 
compared to the second term, Eqn. (102b) for this region will become: 

2 2 

NR = 10 log [(1 - f a ) 2 + 4ll) 2 ( f f ° ) 2 ] 
db yc o yc f 

for f = f Q (104a) 

when f = f Q , Eqn. (104a) becomes: 

NR,, = 10 log (1 - — f a ) 2 = 20 log (1 - — f a) = 

db ° yc o ° yc o 

2 ttp h 

=20 log (1 2- fa) for f = f (104b) 

6 pc o o 
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For the case in Eqn. (104b) where all the terms are positive including 

f >0 and a > 0, one obtains: 
o 


“W' - V < 0 


(105a) 


which is found experimentally also as a negative value. From Eqn. 
(104a) and Eqn. (104b), one finds in the neighborhood of f . f : 


®db (f * V > ®db ( ' ■ 'o' 


(105b) 


«db (f •= V =• “db <£ ’ f o> 


(105c) 


and the curve NR,, (log f) has a local minimum at f = f . Thus, the 
aD o 

value of f Q , the fundamental resonance frequency of the plate, may 

be found from the experimental curve NR^(log f ) . 

Since the fundamental resonance frequency of the plate f can be 

found from the experimental curves NR^(log f), as well as the value of 

NR,, (f = f ) , one is able to calculate the damping factor a of the 
db o 

vibrating plate from Eqn. (104b) as follows: 


10 


, m db (f • V . 

20 = a - is f „) 

yc o 


(106a) 


and since NR,, (f = f ) < 0 and y = ~^r, one has from Eqn. (106a): 

QD o p n 


a = 


yc 

47rf 


[i - 


- NR db (f=f o ) 


■] = 


PC 


2ttp hf 
P o 


[1 - 


-^db^V 


] (106b) 


10 


20 


10 


20 
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where a>0 if NR,, (f = f ) <0. 

db o 

Thus, the values of the fundamental resonance frequency of the 
plate f and the damping factor of the vibrating plate a may be found 
from the experimental curves NR^(log f ) . 

C. The stiffness controlled region : This is the region where 

the frequency of the acoustic wave f is smaller than the fundamental 
resonance frequency of the plate f such that f < f Q . Since |a| << 1, 
the first term in Eqn. (102b) could be neglected and in this region 
Eqn. (102b) becomes: 


2 2 

®db - 10 1( * t(^) 2 (^-^-) 2 ] ■ 
... ,4j f o 2 - f2 | 

= 20 log [— -f ] = 


2irp h If 2 - f 2 1 

" 20 lo * — - f- 


■] = 


r Pp h |f G 2 - f2 

20 log [— ^ 


f — 3 ’ 20 lo * «£> " 


It 2 - f 2 i 

* 20 log j 20 log ( 

P 


(107) 


By taking smaller acoustic wave frequencies f further away from the 

2 2 

fundamental resonance frequency of the plate f , such that f << f Q , 


(107) becomes: 


Eqn. 
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a 2 <• o .2 


^4ir ^9 


NR., = 10 log(^)‘(-^-r- 20 log(^ -*-) 
db yc f pc f ' 


4irf 

= 20 log (-^-) " 20 log f = 


2irp hf 

= 20 log [ £-2-] - 20 log f = 

pc 


= C. - 20 log f 


(108) 


In this region raising the frequency f of the acoustic wave by one octave 

(= doubling the frequency f) will reduce NR^ by 6 decibels. Drawing 

NR, , (log f) on a log-log paper will give a straight line with a slope 
db 

of -20 (decrease of NR^ by 6 db for every octave of f) and a constant 
such that: 

h f o 2 

,1R db (f - 1) 'S' 20 loR < - pc ' > ’ 

2irp h 

= 20 log ( 2-) + 40 log f 

pc o 

= C + 40 log f 
o 

NR db (f=100) = C 1 - 40 db 

Once the fundamental resonance freqeuncy of the plate f has been 

o 

found and the damping factor of the vibrating plate a has been calculated, 

the Noise Reduction coefficient NR,, can be drawn as a function of the 

db 

frequency f of the acoustic wave by using Eqn. (102b) for the region of 
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frequencies where |kdjJ << 1. By drawing the curve NR^Clog f) on a 
log-log paper, one will obtain straight lines in the mass controlled 
region and in the stiffness controlled region. 
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CHAPTER VIII 


FREE VIBRATIONS OF THE FINITE PANEL 

In the present chapter the characteristic resonance frequencies of 
the vibrations of the free clamped finite panel (plate) coupled with the 
acoustic fluid (air) in the rigid duct will be discussed. 

Let the clamped rectangular plate of dimensions "a" and "b" 
be situated in the rigid duct at z = 0 in the x-y plane. The lateral 
dispalcement p(x,y) of the plate will be subject to the boundary conditions 
given in Eqn. (59) as follows: 


n = 0 at x = 0, x * a, y = 0, y = b (109a) 

= 0atx = 0, x = a and ~ = 0 at y = 0, y = b (109b) 


Assuming harmonic time variation e la>t , the partial differential equation 
which governs the lateral displacement of the vibrating panel (plate) 
is given in Eqns. (70) as follows: 


v 4 n 


Y " = D P z 


where: y 


4 p p hw 

* V = — 1 C 


12(l-v 2 )p u 2 

2 

Eh 2 


(1/m 4 ) 


( 110 ) 


p is the external net force per unit area in the positive z-direction 
z 

2 

(N/m ) and the other quantities are defined in Eqn. (70a) . 

Two acoustic waves of mode (m,n) can propagate in the rigid duct of 
dimensions "a" and "b" away from the freely vibrating plate. In 
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accordance with Eqn. (49a), the acoustic wave of mode (m,n) propagating 
in the positive z-direction away from the vibrating plate, will be of the 
forms 


= P cos 
o 


mux 

a 


cos 


nuy i(qz-cut) 
b" 8 


z > 0 (111a) 


where from Eqn. (38b) one has: 


qP 


q o mux niry x(qz-cot) n /......s 

u , = — p . = cos cos — e n z > 0 (111b) 

z+ cop + top a b — 


Similarly, the acoustic wave of mode (m,n) propagating in the negative 
z-direction away from the vibrating plate will be of the form: 


„ mux nuy -x(qzfcot) „ s 

p = -P cos cos — e ' z < 0 (112a) 

- o a b — ' 

-q _ qP o mux nuy -i(qz-hot) „ /.-.o.s 

u = — x p = cos cos — rf- n z < 0 (112b) 

z- <op - cop a b — 


The value of the veolcity of the plate in the lateral direction should 
be equal to the velocity vector of the acoustic wave on either side of 
the plate in accordance with Eqns . (77a) and (78a): 

icon = u z+ (z = 0 + ) = u z _ (z = 0_) (113) 

Because of Eqn. (113), one has the same amplitude P q in both waves in 
Eqn. (Ill) and in Eqn. (112). Substituting Eqn. (111b) and Eqn. (112b) 
in Eqn. (113), one finds the plate lateral displacement n(x,y) in the 
form: 


u 


zp 


_ dn _ 

dt 
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9 A cos 


n^y -lcot 
:os — e 
b 


z * 0 (114) 


iqP 

where A = — - — (m) 

03 p 

The value of the lateral displacement q(x,y) found in Eqn. (114) 
is directly related to the acoustic wave mode (m,n) propagating in the 
rigid duct.. This value of n( x *y) in Eqn. (114) obeys the boundary 
conditions, Eqn. (109b), for a clamped plate at the edges. However, the 
value of p(x,y) at Eqn. (114) does not obey the boundary conditions, Eqn. 
(109a), for a clamped plate at the edges. This requirement that p = 0 at 
the edges of the clamped plate will introduce an interaction between the 
different modes (m,n) given in Eqns . (Ill), (112) and (114) and the funda- 
mental mode (m = 0, n = 0) of the acoustic wave propagating in the rigid 

duct and discussed previously. The value of A in Eqn. (114) for 

m, n 

each mode (m,n) will be related to all the other modes (m,n) and to 
Aq q of the fundamental mode. However, in the present report we are 
interested in finding the values of the characteristic frequencies of 
vibrations and not their amplitudes. We are able to use Eqn. (114) for 
finding these characteristic resonance frequencies 'of vibration, but 
not their amplitudes. 

From Eqn. (111a) and Eqn. (112a) , one obtains at z = 0: 



(p. - P+3 2 =o - -2P 0 


mirx mry -iwt 

cos cos —r^~ e 

a b 


(115) 


Substituting Eqns. (114) and (115 ) in Eqn. (110) , one obtains after 
cancelling the common functions: 


- 2P 

r /Hitt v 2 , ,n*n\ 2 n 2 4 1 o 1 

+ ] ■ Y "d — ’d 


•o 2 

i2cd p 


(116a) 
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, one obtains from Eqn. (116a) : 


2 4 p ha) ' 1 

Substituting y = — ^ and y = — 

P P 


U M)2 + (SI) 2 ] 2 . / 

a d 


4 

i ^ 


(116b) 


For an acoustic wave of mode (m,n) in a rigid duct, one has from Eqn. 
(50b) s 


= \ p ’ - k c 2 =^ k 2 - ( Jf ) 2 - (^) 2 =J £) 2 


(— ) 2 - (- 2 f ) 2 ( 116 c ) 

a b 


Substituting Eqn. (116c) in Eqn. (116b), one has: 


ri 4 4 . [72 , 2 .4 

[k c " Y V k ' k c = lpY 


(117a) 


where k 2 - (— ) 2 + (^-) 2 and k is the cut off wave number for the 
cab c 

acoustic wave mode (m,n) in the rigid duct given in detail in Eqns . (40) 
and (50) . Eqn. (117a) may be rewritten alternatively in the form: 


tV * * V c 


4 4 h /k 2 .k 2 - 4 


yy 


(117b) 


where y = -^rr is the coupling coefficient between the acoustic waves 

v 

and the plate vibrations in the duct, 

4 p o h 2 

Taking k = ~ and y 35 -g- w and substituting in Eqn, (117b), one 
obtains: 


r , 4 ^ 2 f~l XT 2 

[k c “ D “ 1 V k c ‘ ( T } y D 


(118a) 
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By squaring both sides of Eqn. (118a), one obtains: 



2 2 

u) ] [k, 


9 2 « P _ b 9 

2 - —1 = y 2 (- E -) 2 
•c 2 J ^ D ; 

c 


4 

to 


(118b) 


2 

Eqn. (118b) is a cubic algebraic equation in to for the characteristic 
resonance frequencies for mode (m,n) . It has three distinct positive 
roots for the resonance frequencies to of the coupled system of the clamped 
plate and the acoustic waves in the rigid duct. It can be solved for 
numerical values by using the standard cubic equation solution for 
numerical approximation. 

An alternative algebraic solution of Eqn. (118b) for the case 
considered in the present report will be based on the following 
inequalities : 

c 2 » k 2 = [(HI) 2 + (M) 2 ] (119a) 

c a b 



- 

ph 


(119b) 


where in Eqn. (119b) the fundamental mode (m=0, n=0) is excluded. 
Opening the brackets on the left hand side, one of the terms will be 

much larger than the term on the right hand side because of Eqn. (119b) 

2 2 
where k c " » y : 


p h , 

+(^-)' 


, 2 4 , 
k a) + 
c 


» y C 


9 

2 


^ 4 
) oj 


As a result, for the first approximation, the coupling coefficient term 
on the right hand side of Eqn. (118b) could be neglected, and Eqn, (118b) 
becomes for the first approximation: 
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(120) 




D 


2 i 2 r, 2 

“ 1 t k c 



For the first approximation in Eqn. (120) the acoustic resonance fre- 
quency and the plate vibrations resonance frequencies are uncoupled. The 
acoustic resonance freqeuncy can be obtained from Eqn. (120) as 
follows : 

V* - \ - 0i - ck c - c y ( =L)2 + (M)2 (121a) 


where the acoustic resonance frequency for the first approximation 
Eqn. (121a) is equal to the cut-off frequency to = ck c of the acoustic 
wave of the same mode (m,n) in a rigid duct given in Eqn. (50), The 
corresponding plate vibrations resonance frequencies co 2 and can be 
obtained from Eqn. (120) as follows: 



. tiL 

D 



0; o) 2 = id 





+ (S2-) 2 ] (121b) 


where the plate vibrations resonance frequencies w 2 and oj^ for the first 
approximation, Eqn. (121b) , appear as a double root of the equation, 
and is equal to the resonance frequency of the mode (m,n) of the plate 
vibrations given in Eqn. (66). For the first approximation, when the 
coupling effect between the plate and the fluid (air) is neglected, one 
finds the acoustic resonance frequency in Eqn. (121a) as a single 
root of the Eqn. (120), and the plate resonance frequencies oo 2 = in 
Eqn. (121b) as a double root of the Eqn. (120). 
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For the second approximation of cu^, the small coupling term on the 
right hand side of Eqn. (118b) will be included* Substituting the value 
of oj 1 found in Eqn. (121a) in the first and the third terms of Eqn. 
(118b), one will obtain the following equation for the second approxi- 
mation of 


4 , . .2,2., 2 

IK - d (ck c> 1 Ik c 


0) 


1 , 2 , p p h , 2 , . ,4 

— ] - u <-£"> (<*.> 


which may be rewritten in the form: 



2 r p p hc .2 


2 P o hc 2 

fk E 1 

1 c D J 


2 

U 


Dk 


[1 - 


-7 J 


p he 
P 


2 2 

Using Eqn. (119a) where k, « c , one could expand the previous denomi 
nator to give: 


“l 2 2., + 2Dk c 2 

— 5" = U [1 + 2 

C P he 


] 


Rearranging, "one obtains: 


2Dk 

2 2. 2 2 2 , . , c . _ 

a). = c k - c p (1 + o) - 

1 c 


p he 
P 


„ „ 2 2Dk 

c 2 k 2 [1 - -^(1 + ^x - )] 

c k 2 p he 2 

c p 


2 2 

Using Eqn. (119b), where p << k £ , one could expand the previous 
equation to give: 
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(122a) 


“i • ck c [1 ‘ \fi ■ 

k 

c 


2 ■ 2 
Du 


p he 

M P 


2 2 2 

Since from Eqn. *(119) y << k << c , one could neglect the last 

c 

term, as small of the second order of magnitude, for the present appro- 
ximation, and obtain: 

“1 ■ = k c t 1 ■ ck c -liT" 2 <122b) 

k c 

c 


Equation (122b) gives the second approximation for the acoustic resonance 
frequency 0 )^, where the coupling coefficient between the plate and the 
fluid (air) has been included. From Eqn. (122b) , one sees that the 
resonance frequency of the acoustic wave is smaller than the cut-off 
frequency u> c = ck^ of the same mode (m,n) of the acoustic wave propaga- 
tion in the rigid duct: 


U) 

1 c 

a) < ck = c — = oo 
1 c c c 


» in , n * th , n 
f- < f * 

1 c 


(122c) 


Therefore, the acoustic resonance of the plate-air resonance system 
will establish acoustic wave mode (m,n) of frequency in the rigid duct, 
which will not propagate, but will attenuate exponentially as one goes 
further and further away from the free vibrating plate and its coupling 
with the fluid (air). 

For the second approximation of and the small coupling 
term on the right hand side of Eqn. (118b) will be included. Substituting 
the value of = u) » found in Eqn. (121b \ in the second and third 
terms of Eqn. (118b) , one will obtain the following equation for the second 
approximation of and ou: 
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4 , 4 

/ P h 2 - n Uk O P_k ? ^k o 

Pp 


which may be rewritten in the form: 


/ p 1> 2 9 

r k ^ _ _E_ A 
1 c D u 2,3 


2 - 8 
M * 


2 . 6 

!V 


„ Dk 4 Dk 2 

[k c - ~ Sl i) * — rri 


p p hc 


p he 
P 


2 2 

Using Eqn. (119a) where k c << c , one could expand the previous 
denominator to give : 


h 2 


Dk 


[k 


i2 2, 6 r , . c , 

U) 0 ,] = y k [1 + y] 


D 2,3 


p he 
M P 


2 2 

Taking the square root of both sides and expanding, since k c << c , 
one has: 


h 2 


Dk 


—jr~ - =+p k 3 [1 + 

D 2 ’ 3 c 2p he" 

H P 


Rearranging, one obtains: 


p h 2 . , Dk 


2p he 
P 


Dk 


* k c 4[1 4 k- (1 + — 
c c 2p he" 


which will give: 


2,3 


A k e 4[1 1 IT (1 + ~ TT>1 

p c 2p he 
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2 2 

Using the inequality in Eqn. (119b), y << k c , one could expand the 
previous equation to give: 



[1 



(1 + 


Dk 

— ^)] 
2p he 
P 


Taking the positive sign for and the negative sign for one has: 


“2 

w 3 


spL. k 2 ri + I JL 

V P h c [1 + 2 k 


Dk y 

+ Hi 


c 4p he 



-2=-k 2 r, -IJL. 

° P h c 


Dk y 

c 


2 k . , 2 

c 4pphc 


(123a) 

(123b) 


2 2 

Using Eqn. (119a), where k c << c , one could neglect the last term as 
small of the second order of magnitude for the present approximation, 
and obtain: 


^2 



k > + iir ! ■ 

c 



+ 2 ^ k c l 


(124a) 





[1 - i -E 

L 2 k 



(124b) 


Equations (124) give the second approximation for the plate resonance 
frequencies and where the coupling coefficient between the plate 
and the fluid (air) has been included. In this case, one obtains two 
single roots, o) 2 ^ cj^. of Eqn. (118b);, which become a double root, u) 2 = 
only when the coupling coefficient has been neglected. The difference 
in the frequency between the two roots may be found from Eqns . (124) to 
give: 
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(124c) 


a>2 " W3 



The three discrete resonance frequencies oj^, -o^* for the mode 

(m,n) of the coupled system of the clamped plate and the fluid (air) 

in the rigid duct have been given in Eqns. (122b), (124a) and (124b) 

2 2 2 

under the assumption, Eqn. (119), that u << k << c . One could 

c 

rewrite Eqn. (118b) therefore in the form: 


(a ) 2 - oij 2 ) (a ) 2 ~ 1^2 ^ ) (u > 2 


') = 0 


(125) 


2 2 2 

Since ^ ^ , u)^ are the roots of Eqn. (118b). Substituting oj^, 

too from Eqns. (122b), (124a), (124b) respectively in Eqn. (125) 

■ 2 

one obtains Eqn. (118b), if the terms — s- << 1 are neglected, in 

k 1 
c 

accordance with Eqn. (119b) and our present approximation. This checks 
the results found for co^, ^ > W3 under the present approximation. 
Comparing Eqn. (122b) with Eqn. (124), if one had k c << c for the pre- 
sent cases under consideration, one will obtain: 


to 2 << < u> c <<: ( 126 ) 

Thus, the two plate discrete resonance frequencies are below the cut-off 
frequency of the corresponding acoustic wave mode (m,n) in the rigid 
duct and will attenuate exponentially. However, since the plate resonance 
amplitude ^ is coupled with ^ through the boundary conditions 
D = 0 at the edges, as discussed previously, the plate resonance at 
frequencies oj ^ and will produce also the fundamental acoustic mode 
in the duct because of the above interaction, and the fundamental 
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acoustic mode in the duct will propagate at any frequency, regardless 
of the value of the cut-off frequency. 

Let the fundamental mode of the acoustic waves in a rigid duct be 
incident on a panel (plate) situated at z = 0, If one takes the solution 
of the plane wave incident on an infinite plate given in Eqns. (81), one 
sees that the amplitude of the plate vibration is given in Eqns. 

(74a) and (81a), and is a function of time only and not of position on 
the plate. By introducing the duct rigid wall, parallel to the lateral 
direction of the plate, one introduces the boundary conditions, Eqn. (109), 
to the clamped plate, but does not effect the acoustic plane wave, 
which is identical with the fundamental mode. In accordance with Eqn. 
(109a), the clamped plate is not allowed to move at the edges; but 
according to the solution for infinite plate in Eqn. (81a), the plate 

does move at the edge. Since at the edges of the clamped plate n = Q> 

and the clamped plate cannot move at the clamped edges, the plate 
starts to vibrate at higher modes (m,n) in such a way that the total 
net movement of the edges obey q - 0. Thus, it is necessary that the 
clamped plate will vibrate at higher modes in order to compensate for 
the movement of the edges because of the plane acoustic wave, and the 
vibration is such the q = 0 is kept at the edges of the clamped plate. 
Thus, the existence of higher order acoustic modes (m,n) in the duct 

are necessary for the clamped plate case, in order to match the boundary 

conditions q = 0 at the edges of the plate; the higher order modes 
exist even if the incident aocustic wave on the clamped plate is a per- 
fect fundamental acoustic wave, a completely uniform plane wave within 
the duct rigid walls, coming from a perfect source. 
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CHAPTER IX 


THE BERANEK TUBE RECEIVING CHAMBER 

The receiving chamber of the Beranek tube used in the experimental 
set up includes highly nonuniform lossy material of different kinds. 

The total length of the receiving chamber is about 90" inches with 
inside cross-section of 18" x 18". The inside tube contains 18" 
inches of loosely packed foam behind the receiving microphone, with 
additional length of 9" of loosely packed foam in wedge form on both 
sides of the giving microphone. The loose foam on both sides of the 
receiving microphone and behind it absorbs high frequencies in the 
neighborhood of the receiving microphone, and reduces transverse 
standing waves. 

Beheld the foam there are four highly absorbent fiberglass wedges 
with overall length of 60" inches, which include a rectangular base 
24" inches long. The wedges are 9 sq. inch and made from semi-rigid 
fiberglass to help preserve their shape. The density of the fiberglass 

3 

used in the wedges is 3 Ib/ft . The wedges are positioned in the tube 
with their blades alternated to present a maximum absorbent area to the 
sound. The wedges are followed by 8" inches of cotton and 4" inches of 
fiberglass to the back panel. The cotton and fiberglass packed between 
the fiberglass wedges and the back wooden panel help absorb the lower 
frequencies. 

In order to simplify matters from theoretical point of view, it 
will be assumed that the absorption material is uniform and presents an 
average uniform absorption effect to the incident and reflected acoustic 
waves. It will be further assumed that only plane waves propagate in the 
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absorption material. Since the higher order modes attenuate much 
faster than the plane waves in the duct, because of losses at the 
walls of the duct, this assumption should not change much the basic phy- 
sical phenomena. Once the basic absorption phenomena of the incident 
and reflected waves has been analyzed and compared with experiment, 
additional corrections can be made by taking into account the nonunifor- 
mity of the absorption material and the higher order modes of propagation. 

Let a uniform plane acoustic wave be incident in the positive z- 
direction in a lossy medium: 

-<xz i(kz-oot) (ik-<*)z-ioot /, 0 7 ^ 

p ± = p. e d = P I e (127a) 

where k » oo/c and « represents the attenuation factor in the lossy 
medium. Substituting Eqn. (127a) in Eqn. (20b), one obtains: 


ik-°c — ocz i(kz-cot) l+i«/k +(ik-«) z-icot 

u . = p T - e e = p_ e 

Zl r I loop I pc 


(127b) 


Let a uniform plane acoustic wave be reflected in the negative z-direc- 
tion in a lossy medium: 


p r = p R e • 


-i(kz-hi)t) _ -(ik-*)z-iojt 


= P R e 


(128a) 


Substituting Eqn. (128a) in Eqn. (20b) , one obtains: 


(ik-°0 a z -i(kz+u>t) 14-i^/k -(ik-oc) z-ioot ^ noi X 

u - — p-, : e e - -p^, e (,1/ob; 

zr R loop R pc 
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Taking the back wooden panel to be situated at z = L, one has there 

u ■ 0 and from Eqns. (127b) and (128b) one obtains, after cancelling 
n 

the common factors: 


u = 0 = u . 
n zi 


, (ik-«)L -(ik-*)L ,, OQ . 

+ u zr I P I e " P R e (129a) 

z=L z=L 


From Eqn. (129a) one has 


P R = Pi e 


2 (ik-®)L 


(129b) 


which is identical with Eqn. (28b) for « = 0. Assuming in the present 
case that the attenuation factor is small, °=/k << 1, one is able to 
continue the analysis similarly to Eqns. (29) to (35). The effect of 
the attenuation factor * ^ 0 is to reduce the value of the resonance 
frequencies calculated for « = 0. 

Rewriting Eqn. (35), one has: 

P max = z 0 ) = P I |cos[k( L -z 0 ) * (to) ] | (130) 

For the case d> (uj) = 0, one requires for p = 0 at z = z that: 
r ^ r max o 


k(L - z Q ) = f (L 


Z 0 } = 


2TTf 


a - z 0 ) = 


2ir 


(L 


, _ (2n + l)j (131a) 
z 0 ) - 2 


where n = 1,2,3 . . . From Eqn. (131a) , one obtains the characteristic 
frequencies for minimum value of the pressure p at the receiver microphone 
situated at z = Zq? 
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f n“ WqT 1 for ♦ <»> = 0 (131b) 

For the case c|> (to) = ~7r t one requires in Eqn. (130a) for p = 0 at 

/ max 

z = Zq that: 

k(L - z 0 ) + f = ~ (L - z q ) + | = (2n + 1) | (132a) 

where n = 1,2,3 . . . From Eqn. (132a), one obtains similarly: 

for -Kw)* 1 ? (132b) 

n Ll — Zq £ £ 

In our experimental set up the distance of the wooden back panel from 
the test aluminum panel is L = 107.5" and the distance of the receiver 
microphone is Zq = 8". Therefore, one has in our experiment: 

L - z = 107.5" - 8.0" =99.5" = 2.527 m (133a) 

o 

Taking the acoustic velocity c for T = 21°C = 70°F from Eqn. (6) 
to be : 


c = 343.8 m/sec (133b) 


Using Eqn. (133) in Eqns. (131b) and (132b), one obtains: 

f = 68.0(n + h (1/sec) for <f>(u) =0 (133c) 

n z 

f = 68.0 n (1/sec) for $ (id) = ~ (133d) 

n Z 

for n = 1,2,3 . . . 
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From the experimental results it is found that Eqn. (133d) represents the 
results much better than Eqn. (133c). Thus, a phase shift of cf> (oo ) = -90° 
is introduced between the incident plane wave and the reflected plane 
wave by the wooden back panel in accordance with Eqn. (33). As indicated 
earlier, the attenuation factor a ^ 0 in the lossy medium will introduce 
a shift of the actual resonance frequencies found experimentally to 
lower frequencies than calculated in accordance with Eqn. (133d), 
which is a well established phenomena. 

Let us now calculate the resonance frequency of the first even 
vibration mode (m - 2, n = 2) of the wooden back panel in accordance 
with Eqn. (55c) and Eqn. (66d), which are reproduced here for reference: 


f wood 

m,n 


= J [(f ) 2 + (§) 2 ] 



Eh' 


12(1 - v ) 


(134a) 


Taking for our experimental set up requirements a = b and m = n, one 
obtains from Eqn. (134a): 


^wood 
m, m 



h 



_E 

(1 - v 2 ) 


(134b) 


In the present experimental set up the wooden back panel has the follow- 
ing parameters: 



18" = 0.4572m; h = 0.75" = 0.019m 


6,89 lb -, = 49.0 lb/ft 3 = 785.0 Kg/m 3 
243 inch 


(135a) 

(135b) 
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'Since the Young's modulus of elasticity E is not available for the 
wooden back panel, it will be assumed to be identical with poplar 
wood as follows: 

E = 1.0 x 10 10 (N/m 2 ) 1 - v 2 = 0.9 (135c) 

Substituting Eqns. (135) in Eqn. (134b), one obtains the resonance 
frequency of the first even vibrational mode m = 2 of the wooden 
back panel: 


.wood 

2,2 


= 1,240 1/sec 


(135-d) 


The experimental curve in the present paper shows a very sharp 
resonance frequency at f - 1,260 Hz which could be identified with the 
above vibrational mode m = 2 in Eqn. (135d) of the wooden back panel. 

The above resonance represents interaction between the incident acoustic 
and reflected acoustic waves with the vibrational modes of the 
aluminum panel in front, the wooden back panel and the receiver micro- 
phone. A more detailed analysis of the wooden back panel will be 
given in the next chapter. 

When two resonance systems interact with each other the resonance 
frequency oo of the combined system could be given similarly to Eqn. (125) 

(O) 2 - I^ 2 ) (U3 2 - 0) 2 2 ) = y 0 4 (136a) 

where is the resonance frequency of the first system, is the 
resonance frequency of the second system (w > w^) and represents the 
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coupling factor between the systems . When there is no coupling 
Uq = 0, the resonance frequencies are: 

Vig = 0; to = io_ = co^ and to = co + = > 10 ^ (136b) 

In case of coupling between the systems 0, Eqn. (136a) can be 
rewritten in the following form: 

m 4 - Go/ + co 1 2 )(o 2 + (co 2 2 ( o 1 2 - v Q 4 ) = 0 (137a) 

2 

Solving the quadratic Eqn. (137a) for <0 for ^ 0, one has: 

to/ = \ [(co 2 2 + CO, 2 ) + ^/(co 2 2 - to/) 2 + 4y 0 4 ] > < 0 / (137b) 

to/ = \ [(co 2 2 + co 1 2 ) - 0 / - to/) 2 + 4y 0 4 ] < to/ (137c) 

where Eqns. (137b) and (137c) represent the natural resonance fre- 
quencies of the coupled systems, when the oscillators are no. longer 
isolated. Thus, the coupling between the systems causes the higher 
resonance frequency io + to shift higher, the lower resonance frequency 
to to shift lower. When there are several resonance systems interact- 
ing with each other, it is not surprising sometimes to find that the 
experimental resonance frequencies are shifted as compared to the 
theoretical resonance frequencies calculated for each system by 
itself, without taking into account the coupling factor Pq f 0. 
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CHAPTER X 


EXPERIMENTAL RESULTS 

In the present chapter we will apply the theory which has been 
given in this report to analyze some of the outstanding apsects of 
the experimental results given in Figure 1. The aluminum panel which 
gives the noise reduction in Figure 1 has the following dimensions and 
properties : 

Material: Alclad 2024T3 Aluminum 

Thickness: h = 0.025" = 0.635 mm 

Width x Axis: a = 18" = 0.4572 m 
Width y Axis: b = 18" = 0.4572 m 

Some mechanical properties of aluminum are as follows: 

3 3 

Density = p p = 2.7 x 10 Kg/m 

10 2 

Young modulus of elasticity = E - 7.0 x 10 N/m 
Poisson's ratio = v =0.3 

1 - v 2 = 0.91 

The experiments were done in Lawrence, Kansas at about 1,000 ft above 
sea level under the following conditions of the air: 

Temperature - T « 21°C = 70°F = 294.15°K 

2 

Pressure = P = 0.97735 atmosphere = 97,735 N/m 
Density = p = 1.1575 Kg/rn^ 

Velocity of Acoustic Waves = c = 343.8 m/sec 
The analysis of the experimental results in Figure 1 will be 
divided to several parts according to the particular aspect of the phy 
sical phenomena to be discussed. 
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A. Plane Wave on the Finite Panel 


Since for the transmitter microphone = 1 cm one has |kd^| « 1, 
the general characteristics of the Noise Reduction curve of Figure 1 
could be found by using Eqn. (ll>2b) which is rewritten here: 


f 2 f 2 f 2 

, r / 1 4ir 0 a ,2 , /4 tt n 2 , ■ 0 ,2, 

NR db = 10 log[(1 - y c~~ ) + W ( f } ] 


(138) 


The coupling parameter p between the aluminum panel and the air is: 


_ 2p _ 2 x 1.1575 _ - , , 

y p h 2.7 x 0.635 1,350 (1/ ^ 

P 


(139) 


From Figure 1 one sees that the fundamental resonance frequency of the 
plate is: 


fg = 63 Hz (140a) 

at which point the Noise Reduction is: 


NR db (f = f Q ) = -4.5 db 


(140b) 


Using the experimental values of Eqn. (140) in Eqn. (106b), one is 
able to calculate the damping factor a by using Eqn. (139) and 
c = 343.8 m/sec: 
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a - Jf- [1 - 

4irf„ 




1.350 x 343.8 
4 it x 63 


U - -frl = 


10 


20 


10 


*tA 

20 


= 0.58626 [1 = 0.58626 [1 - -- ™o" Q 3 = 0-2370 


10 


0.225' 


1.6788 


Substituting Eqns. (139) and (140) in Eqn. (138), one obtains; 


„„ , r/, 25.468 v 2 , ,^-4,f^ - 63^2, 

NR db = 10 log [(1 j ) + 7.33 x 10 ( j ) ] 


where NR^ is given as a function of the frequency for Figure 1. 
gives the results for some of the frequencies calculated from Eqn 


TABLE B 

NOISE REDUCTION VALUES OF EQN. (141) 


f 

M db 

20 Hz 

+13.7 db 

30 Hz 

8.9 db 

40 Hz 

+4.3 db 

50 Hz 

- 0.6 db 

63 Hz 

- 4.5 db 

70 Hz 

- 2.7 db 

80 Hz 

+ 0.6 db 

100 Hz 

5.1 db 

200 Hz 

13.9 db 

400 Hz 

20i.5 db 

1000 Hz 

28.7 db 

5000 Hz 

42.6 db 


(140c) 

(141) 

Table B 
. (141). 







Figure 1 gives the experimental Noise Reduction curve for the 
aluminum panel discussed here. The thin line gives the experimental 
ratio of pressures in decibels between the receiver microphone located 
8" behind the aluminum panel to the transmitter microphone located 
1 cm in front of the aluminum panel, as a function of frequency. The 
thin straight line in Figure 1 represents a straight line least square 
error average of the experimental curve. The heavy line represents 
the theoretical curve of the present paper, with the values given in 
Table B and calculated in accordance with Eqn. (141). As shown above, 
the theoretical curve is in accordance with the experimental 
fundamental frequency of resonance f^, and the experimental Noise 
Reduction value at this frequency, so that the theoretical curve is 
passing through this point f = f^ on the experimental curve. 

The difference between the experimental curve (thin wiggly line) 
and the theoretical curve (thick black line) represents two major 
effects : 

1. In the higher frequency region this difference represents the 
effect of the higher order modes and the different resonance 
phenomena . 

2. In the lower frequency region this difference represents the 
effect of the near zone fields next to the panel. 

B. The Fundamental Resonance Frequency (f^ = 63 Hz) 

The fundamental experimental resonance frequency f^ =63 Hz in 
Figure 1 and in Eqn . (140a) is the result of interaction between two 
major resonance phenomena: 
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1. The standing plane acoustic wave which propagates along 

the length of the Beranek tube and is reflected by the wooden 
back panel. This basic resonance frequency is given in Eqn. 
(133d) f^ avit y = 68 Hz. However, this resonance frequency 
is shifted downward because of its interaction with the vibra- 
tion resonance frequency of the wooden back panel. 

2. The first even resonance mode of the aluminum panel 

vibrations, which excite acoustic waves at = 58.8 Hz 

z> 9 z 

(see later). Although this acoustic wave mode does not pro- 
pagate in the high pass filter duct, because of its large 
wavelength it has a definite effect on the receiver micro- 
phone as a near zone acoustic wave field excited by the first 
even mode of the aluminum panel vibrations. 

As a result of the interaction between the above two basic resonance 
effects, one obtains the fundamental experimental resonance frequency at 
fg = 63 Hz. 

The evidence for the above interaction can be found in other experi- 
mental data: 

a. From Eqn. (132b) it is found that f^ avlt ^ ± s inversely 
proportional to the total length L of the Beranek tube 
since |zq| << L. Experimental data has shown that by 
decreasing the length of the Beranek tube the fundamental 
resonance frequency f q increases . The same will happen 
even if there is no wooden back panel and the duct end 
is open to infinite space. The open duct end causes 
reflection of the acoustic wave as much because of the 
sudden change of impedance from the duct to infinite space. 
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b. From Eqn. (134b) one finds that the aluminum panel 
resonance frequency of the first even mode f^^ te 
is directly proportional to the thickness of the plate h. 
However, experimental data has shown that the fundamental 
resonance frequency f Q does not change by much for 
changing thickness of the aluminum panel. This is because 
the resonance effect (a) of the standing plane 
acoustic wave along the length of the Beranek tube does 
not change, and it is very dominant in the interaction 
between the two resonance phenomena. 

C. The Cavity Resonance 

The resonance effect of the standing plane acoustic wave along the 
length of the Beranek tube will be called the cavity resonance effect 
and will be discussed in the present section. The cavity resonance 
effect is described by Eqn. (133d) as follows: 

f cavity „ 68 g (s = 1 , 2,3 . . .) (142) 

In the following Table C the theoretical numerical values of Eqn. (142) 
as well as the corresponding resonance frequency experimental 
values from Figure 1 will be given. Theoretical numerical values from 
Eqn. (142) with no corresponding experimental values from Figure 1 will 
be omitted from Table C. 
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Note 













TABLE C (continued) 


s 

jcavuy „ 68 
s 

Experimental 

*30 

2,040 

*2,070 

32 

2,176 

2,200 

38 

2,584 

2,600 

43 

2,924 

2,900 

47 

3,196 

3,180 

*52 

3,536 

*3,500 

55 

| 3,740 

3,760 

*59 

! 

| 4,012 

i 

*4,000 

60 

i 

4,080 

4,070 

65 

4,420 

4,400 

*66 

4,488 

*4,480 

72 

4,896 

4,860 


Note: * Denotes major experimental upward spikes. 


Some of the experimental values listed in Table C are only a 

sharp break in the experimental curve. Most of the experimental values 

indicated in Table C are upward spikes in the experimental curve given 

in Figure 1- The resonance frequency f^ avlt y in Eqns. (142) and (133d) 

are derived from Eqn. (130) under the condition that p (z = z n ) =0, 

n r max 0 

which means that the pressure at the receiver microphone is zero. 
Therefore, at the resonance frequencies listed in Table C, the resonance 
spikes in Figure 1 are upward, representing a strong noise reduction, 




2 











since the transmitted acoustic wave at the receiver microphone is taken 
to be zero at these frequencies. However, because of the near zone 
fields of the higher order modes, the pressure ,at the receiver micro- 
phone is not quite zero, and therefore the resonance upward spikes of the 
experimental curve in Figure 1 are finite. 

The agreement between the theoretical cavity resonance frequencies 

and the experimental resonance frequencies in Table C is relatively 

good since no interaction effect has been taken into account in the 

theory, such effects usually shift the calculated resonance frequency. 

The interaction effects include the near zone fields of the higher order 

modes excited by the panel, and the interaction effects with the wooden 

back panel and the wooden walls of the duct, the effects of which have 

been neglected by assuming that the wooden walls are infinitely massive 

and do not vibrate. Particularly large sharp upward spikes of this 

phenomena appear at the f cav: *- t y for s = 3, 9, 15, 16, 19, 25, 30, 52, 59, 

s 

66 s which are denoted by a star at Table C. 

D . The Acoustic Resonance 

It has been shown in Chapter VIII that the free vibrations of the 
panel will establish acoustic wave resonance inodes (m,n) of circular 
frequency in the rigid duct, the value of which is given in Eqn. 

(122a) and Eqn. (122b) . This acoustic wave resonance inodes will not 
propagate in the duct, but will attenuate exponentially as one goes 
further and further away from the free vibrating plate coupled with the 
air system, because the resonance frequencies of these modes are 
slightly below the cut-off frequency of the corresponding higher order 
modes in the duct, as shown in Eqn. (122c) . Thus, at these acoustic 
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- resonance frequencies of the free vibrating plate, the acoustic pressure 
from these modes will be higher air- the transmitter microphone than at 
the receiver microphone, since the transmitter microphone is much 
closer to the aluminum panel than the receiver microphone. As a 
result, at the acoustic resonance frequencies one will get upward 
spikes in the experimental curve representing resonance "noise reduction." 

Equation (121a) gives the acoustic resonance frequency in the 
following form for a = b: 


acoustic 

1 


w, ck i- - 

1 c c /.mir.2 . .mr.2 

- IT - V ( T> + ( T> 




(143a) 


Substituting c = 343.8 m/sec and a = 18" = 0.4572 m, one obtains from 
Eqn. (143a) for the acoustic resonance frequencies: 


^acoustic 

m,n 


375.98 




376.0 


2 2 
m + n 


(143b) 


In the following Table D the theoretical numerical values of Eqn. (143) 
are arranged in order of increasing acoustic resonance frequencies f^ C ° 
The corresponding resonance frequency experimental data of the upward 
spikes are given from Figure 1. 
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TABLE D 


ACOUSTIC RESONANCE FREQUENCIES 


m - n 

2 . 2 
m + n 

^acoustic 

m,n 

Experimental 

*1-0 

1 

376.0 

*369 

1-1 

2 

531.7 

517 

**2 - 0 

4 

752.0 

**742 

2-1 

5 

840.7 

798 

*2-2 

8 

1,063.4 

*1,080 

*3-0 ’ 

9 

1,127.9 

*1,080 

3-1 

10 

1,189.0 

1,220 

*3-2 

13 

1,355.6 

*1,300 

4-0 

16 

1,503.9 

1,470 

4-1 

17 

1,550.2 


3-3 

18 

1,595.1 

1,600 

*4 - 2 

20 

1,681.4 

*1,720 

4-3 

25 

1,879.9 

1,830 

5-0 

25 

1,879.9 

1,830 

5-1 

26 

1,917.1 

1,940 

*5-2 

29 

2,024.7 

*2,070 

4-4 

32 

2,126.9 


5-3 

34 

2,192.3 

2,200 

6-0 

36 

2,255.9 


6 - 1 

37 

2,287.0 



Note: * Denotes major experimental upward spikes 


^acoustic 

m,n 


376 


V' 


' 2 , 2 
m + n 
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m - n 

2 , 2 
m + n 

^acoustic 

m,n 

Experimental 

! 7-7 

98 

3*722.0 

■ ! 1 


8-6 

100 

3,759.8 

3,760 

10 - 0 

100 

3,759.8 

3,760 

10-1 

101 

3,778.6 

: 

3,760 

10-2 

104 

3,834.2 


'9-5 

106 

3,871.1 


10 - 3 

109 

3,925.2 


*8-7 

113 

3,996.7 

*4,000 

*10-4 

116 

4,049.3 

*4,000 

9-6 

117 

4,067.0 

4,070 

11-0 

121 

4,135.8 

4,070 

11-1 

122 

4,152.7 


10 - 5 

125 

4,203.5 


11-2 

125 

4,203.5 



8-8 128 


9 - 7 

130 

11 - 3 

130 

10-6 

■ ■ 

136 

11-4 

137 

*12 - 0 

144 

*9 - 8 

145 

*12 - 1 

145 

11-5 

146 

12 - 2 

148 

10-7 

149 


4, 253.8 


4, 286.9 
4,286.9 

4.384.7 

4.400.8 

4.511.8 

4,400 

4,400 

*4,480 

4,527.6 

*4,480 

4,527.6 

*4,480 

4,543.0 

j 

4,574.2 


4,589.6 


















TABLE D (continued) 

















Particular large sharp upward spikes of this phenomena appear at 

acoustic 

f m,n for (m ’ n) = (1 ’ 0) ’ (2 ’ 0) » (2 » 2) » 0,2), (4,2), (5,2), 

(7,6), (9,2), (8,5), (8,7), (10,4), (12,0), (9,8), (12,1) which are 
denoted by one star at Table D. Particular prominent large sharp 
upward spikes appear at (2,0) and at (7,6), (9,2), (8,5) which are 
denoted by two stars at Table D. The first large sharp spike at 
(2,0) of the acoustic resonance frequency is explained as excited 
by the basic even particularly large vibrations modes (2,0) and (0,2) 
of the aluminum panel. The other particularly large sharp upward spike 
of the acoustic resonance is explained as the interaction of the 
acoustic modes (7,6), (9,2), (8,5) with the close by cavity resonance 
mode s = 52. 

E. The Plate Resonance 

It has been shown in Chapter VIII that the free vibrations of the 
panel will establish plate acoustic wave resonance modes (m,n) of cir- 
cular frequency ar *d the rigid duct, the value of which is 

given in Eqn. (123) and Eqn. (124). These acoustic wave plate resonance 
modes have a frequency which is much below cut-off of the corresponding 
higher order modes in the duct, and will not propagate at all. Their 
effect is primarily a near zone effect near the plate. These plate 
resonance frequencies will also interact with the plane acoustic wave 
excited by the plate, and this plane acoustic wave could propagate in 
the duct at all frequencies since its cut-off frequency in the rigid 
duct is zero. By these two aspects, the plate resonance modes will 
affect the microphones. The experimental results in Figure 1 show that 
this plane acoustic wave, which results from the interaction with the plate 


resonance modes, is added iti the transmitter microphone out of phase to 
the incident acoustic wave, and is added in the receiver microphone 
in phase to the incident acoustic wave. As a result, the pressure 
in the receiver microphone will be higher than usual at plate resonance 
frequencies, and the pressure at the transmitter microphone will be 
lower than usual at plate resonance frequencies. This will result in 
less "noise reduction" in the experimental data at plate resonance fre- 
quencies, and thus will result in downward spikes at plate resonance 
frequencies in the experimental data given in Figure 1. 

Equation (121b) gives the plate resonance frequencies in the 
following form for a = b: 


- _ - 1 / D . 2 1 / D r .ittiFv 2, ,mr N 2, 

f 2 - f 3 m U V; P'7 k <= 2^h [< T> +( T> 1 


_ TT / D , 2 , 2. 

+n) 


(144a) 


Substituting the value of D = 


Eh' 


12(1 - v ) 


in Eqn. (144a), one obtains: 


-plate _ -plate _ trh / 

-2 " f 3 -' 2a 2V- 


12p p (1 - v ) 


2 2 
— (m + n ) 


(144b) 


Using the numerical values for the present case from the beginning of 
this chapter, one obtains for the plate resonance frequencies: 

f plate = y 353 (tn 2 + (144c) 

m,n 

In the following Table E the theoretical numerical values of Eqn. (144) 
are arranged in order of increasing plate resonance frequencies 


110 


^2m n = ^3 m n * corres P onc ^^ n § resonance frequency experimental 

daita of the downward spikes are given from Figure 1. The resonance 
frequency where particularly large downward spikes are found in Figure 1 
are denoted by a star. The even modes of the plate vibrations 
are found most prominent in (2,2), (4,4), (6,4), (8,8) of the first 
few modes. The plate resonance mode (6,6) downward spike is 
cancelled by the acoustic resonance upward spike of mode (1,1). The 
plate resonance mode (10,10) downward spike is cancelled by the acoustic 
i resonance upward spike of mode (4,0). 


TABLE E 

PLATE RESONANCE FREQUENCIES 


m - n 

2 2 
m + n 

f plate 

m,n 

Experimental 

J 

2-0 

4 

29.4 


*2 - 2 

8 

58.8 

*63 

4-0 

16 

117.6 

102 

CM 

1 

20 

147.1 

134 

*4 - 4 

32 

235.3 

*219 

6-0 

36 

264.7 

247 

6-2 

40 

• . ' 

294.1 

290 

*6-4 

. 

52 

382.4 

*383 

8-0 

64 

470.6 

454 

8 - 2 

68 

500.0 



Note: * Denotes major experimental downward spikes 

fPlate = ?<353 ( 2 + n 2 

m,n 
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TABLE E (continued) 


! 

m - n 

" 

2.2 
m + n 


16-6 

292 


| 14 - 10 

296 


16-8 

320 


18-0 

324 


18-2 

328 


*14 - 12 

340 


*18 - 4 

340 


16 - 10 

356 


*18 - 6 

360 


18-8 

388 


14 - 14 

392 


16-12 ■ 

400 


20-0 

400 


20-2 

404 


*20 - 4 

416 


18 - 10 

424 


20-6 

436 


*16 - 14 

452 


20-8 

464 


18 - 12 

468 


22-0 

484 


22 - 2 

488 


20 - 10 

500 


22-4 

500 


16 - 16 

512 
















TABLE E (continued) 



m - n 

2,2 
m + n 

jplaite 

ra,n 

Experimental 

18 - 14 

520 

3,823.6 

3,700 - 3,800 

22-6 

520 

3,823.6 

3,700 - 3,800 

*20 - 12 

544 

4,000.0 

*3,920 

22-8 

548 

4,029.4 

4,060 

24-0 

576 

4,235.3 


18 - 16 

580 

4,264.7 


24-2 

i 

580 

4,164.7 


*22 - 10 

584 

4,294.2 

*4,320 

*24 - 4 

592 

4,353.0 

*4,320 

20 - 14 

596 

4,382.4 

. 

*24 - 6 

612 

4,500.0 

*4,510-4,600 

*22 - 12 

628 

4,617.7 

*4,510 - 4,600 

24-8 

640 

4,705.9 


18 - 18 

648 

4,764.7 


20 - 16 

656 

4,823.6 


24 - 10 

676 

4,970.6 

4,930 

26-0 

676 

4,970.6 

i 

22 - 14 

680 

5,000.0 


26 - 2 

680 

5,000.0 


26 - 4 

692 

5,088.3 

; 

26-6 

712 

5,235.3 


24 - 12 

720 

5,294.2 


20 - 18 

724 

5, 323.6 
















It has been shown in Eqns. (123) and (124) that the plate resonance 
frequencies are split ^2 ^ ^3 a STna -^ va l ue to two components, 
similar to the split of the spectrum of the atom by applying magnetic 
fields. As may be seen from Eqn. (124c) the difference between the 

frequencies P ro P ort ional to k c and becomes larger for high mode 

numbers (m,n). Thus, we see in Figure 1 that as a result, the downward 
resonance spikes resulting from the plate resonance frequencies are 
duller in most cases than the upward spikes, becoming duller for higher 
mode numbers and higher resonance frequencies. Because of the logarith- 
mic scale, narrowing the distinction at higher frequencies, some 
of the higher frequencies downward spikes will look much duller at 
higher frequencies under a regular scale of frequencies. 

It should be pointed out that because of the interaction of the 
above three major effects, the cavity resonance, the acoustic resonance 
and the plate resonance, sometimes the different resonance effects 
add out of phase and cancel each other, leaving no spikes in Figure 1, 
but only a break in the resonance curve, if at all. Thus, it is 
sometimes impossible to detect a particular frequency in the experi- 
mental results of Figure 1, because it has been cancelled by other 
resonance phenomena. Sometimes these different resonance phenomena 
interact with each other to produce a resonance in the experimental 
curve in Figure 1 at frequencies which can not be identified with any 
one of the single phenomena discussed above. 
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F . The Wooden Back Panel Resonance 


Although the wooden back panel is relatively massive, it can act 
as a resonator under the influence of the plane acoustic wave incident 
upon it. Thus, the wooden back panel excites reflected plane 
acoustic waves at frequencies of its own wooden "plate" resonance, 
similarly to the phenomena discussed in the previous section. Rewrit- 
ing Eqn. (144b): 


^wood 

m,n 


= _i*L . / E 

2a 2 Vi2 Pp (l - v 2 ) 


, 2 . 2 , 

(m + n ) 


(145a) 


and substituting the numerical values for the wooden back panel given 
in Eqn. (135), one obtains: 


^wood 

m,n 


155 (m 2 + n 2 ) 


(145b) 


In the following Table F the theoretical numerical values of Eqn. 
(145) are arranged in order of increasing wooden back panel "plate" 
resonance frequencies f^°^. The corresponding resonance frequency 
experimental data of the downward spikes are given from Figure 1. 
Particularly large downward spikes are denoted by a star. 
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TABLE F 


WOOD RESONANCE FREQUENCIES 


1 

j m - n 

2 , 2 
m + n 

f 

m,n 

Experimental f^ n 


4 

620 

*630 

**2 - 2 

8 

1,240 

**1,260 

**4 - 0 

16 

2,480 

**2,420 

*4-2 

20 

3,100 

*3,010 

■HUM 

32 

4,960 

4,930 


36 

5,580 



Note; * Denotes major experimental downward spikes 

f = 155 (m 2 + n 2 ) 
m,n 

wood 

Of particular interest are the large downward spikes of f- „ 

2. yi. 

and f^°°^. The resonance frequency f^ 0 ^ = 1»240 Hz does not have any 
plate resonance frequency close by to interact strongly with it, and 
therefore if. is a sharp downward spike. The resonance frequency 
^4°0^ = 2,480 Hz interacts strongly with the aluminum plate resonance 
frequencies f^ a ^® = = 2,500 Hz and therefore the downward spike 

in this case is a dull one. 
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CHAPTER XI 


SUMMARY 

In the present report the theoretical background has been given and 
the theory has been developed for acoustic plane waves normally inci- 
dent on a clamped panel in a rectangular duct. The coupling theory 
between the elastic vibrations of the panel (plate) and the acoustic waves 
propagation in infinite space and in the duct has been considered in 
detail. The coupling theory developed in this report is based on the 
theory of acoustic waves propagation and the dynamic theory of elasti- 
city and the plate vibrations. The partial differential equation which 
governs the vibrations of the panel (plate) has been modified by adding 
to it stiffness (spring) forces and damping forces, and the fundamental 
resonance frequency £ and the attenuation factor a are discussed in 
detail. 

The noise reduction expression based on the present theory has 
been found to agree well with the corresponding experimental data of a 
sample aluminum panel in the mass controlled region (f > f Q ), the 
damping controlled region (f ~ f' ) and the stiffness controlled 
region (f < f Q )* All the frequency positions of the upward and downward 
resonance spikes in the sample experimental data are identified theo- 
retically as resulting from four cross-interacting major resonance 
phenomena: The cavity resonance, the acoustic resonance, the plate 

resonance, the wooden back panel resonance, and detailed tables are 
given for the values of these resonance frequencies for each case. 

In Chapter I the basic general theoretical considerations are 
introduced. 
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In Chapter II some basic concepts of the theory of propagation of 
acoustic waves and the partial differential equations which govern 
the propagation of acoustic waves in three dimensions are given. 

In Chapter III the acoustic uniform plane wave which varies 
harmonically with time is given using complex phasor representation. 

The incident, ref lected and standing acoustic waves as well as the re- 
flection coefficient are discussed. 

In Chapter IV the porpagation of acoustic waves in rectangular 
ducts is discussed in detail and some of the basic concepts are intro- 
duced. The fundamental acoustic mode, which behaves like a plane 
acoustic wave in infinite space, propagates in the duct at all frequencies 
Higher order acoustic modes attenuate exponentially at frequencies 
below the corresponding cut off frequencies, and propagate at frequencies 
above the cut off frequencies. A list of the cut off frequencies for 
a duct with a square cross-section is given. 

In Chapter V some basic concepts of the vibration theory of elastic 
plates, as based on the theory of elasticity, and the partial differen- 
tial equation which governs such vibrations, are discussed. The 
different possible boundary conditions, and the characteristic fre- 
quencies of vibration are given. 

In Chapter VI the boundary value problem of acoustic plane wave 
normally incident on an infinite vibrating plate is solved, and the 
transmission and reflection coefficients of the corresponding acoustic 
waves are found. The transmission loss coefficient, the noise reduction 
coefficient, and the insertion loss coefficient are defined and are 
related to each other. The solution gives the "normal incidence mass law. 
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This theoretical result agrees with the experimental results in the mass 
controlled region of high frequencies, but does not explain the experi- 
mental results in the damping contolled region (near the fundamental 
resonance frequency) or in the stiffness controlled region (frequencies 
below the fundamental resonance frequency) . 

In Chapter VII the boundary value problem of an acoustic plane wave 
normally incident on a finite clamped panel is solved. The partial 
differential equation which governs the vibrations of the plate (panel) 
is modified by adding to it stiffness (spring) forces and damping 
forces, and the fundamental resonance frequency is defined. The trans- 
mission and reflection coefficients of the corresponding acoustic waves 
are found. The transmission loss coefficient, the noise reduction 
coefficient and insertion loss coefficients are derived. This theoretical 
result agrees with the experimental results in the mass controlled 
region (frequencies above the fundamental resonance frequency) as well 
as in the damping controlled region (near the fundamental resonance 
frequency) and in the stiffness controlled region (frequencies below 
the fundamental resonance frequency). A formula is given for the cal- 
culation of the damping factor from the experimental curve value at the 
fundamental resonance frequency. 

In Chapter VIII the resonance frequencies excited by the free 

vibrations of the finite plate (panel) and its coupling with the air 

2 

surrounding it is discussed. A cubic equation for (27rf) , where f are 
the resonance frequencies for each mode (m,n) of vibration, is solved. 

One solution is found to be the acoustic resonance frequency f^, 

which is slightly smaller than the cut-off frequency of the corresponding 
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mode in the rigid duct, and therefore it will not propagate, but will 
attenuate exponentially. The other two solutions are found to be the 
plate resonance frequencies f ^ ~ f^, which differ in magnitude 
slightly, and are very much smaller than the cut-off frequency of 
the corresponding mode in the rigid duct. These two close together 
plate resonance frequencies will excite higher order acoustic, modes 
which will effect the experimental results only in the near zone, and 
will also interact with the plane acoustic wave fundamental mode, 
which can propagate in the rigid duct at any frequency. 

In Chapter IX the propagation of the plane acoustic wave in the 
receiving chamber of the Beranek tube is discussed. It is shown that 
the reflection of the acoustic plane wave from the wooden back panel will 
cause standing waves and the cavity resonance phenomena. The vibrations 
of the wooden back panel as a vibrating plate is discussed. The 
effect of the coupling between two resonance systems on the resonance 
frequency of each system by itself is considered. 

In Chapter X the experimental results are discussed in detail in 
view of the theory presented in this report. Six major theoretical 
apsects of the experimental results are presented: 

A. Plane Wave on the Finite Panel 

The theoretical noise reduction curve is calculated and drawn in 
heavy line in Figure 1 and follows closely the average experimental 
data in the mass controlled region, the damping controlled region and 
the stiffness controlled region. The difference between the wiggly experl 
mental curve and the heavy line theoretical curve represents the effects 
of the resonance phenomena and the higher order modes. 
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B. The Fundamental Resonance Frequency (f Q = 63 Hz) 

The fundamental resonance frequency f Q = 63 Hz is the result of the 

cavi tv 

interaction between the basic cavity resonance frequency = 68 Hz 

of the standing plane wave along the length of the Beranek tube, and 
the first even resonance mode of the aluminum panel vibrations, 
which excite the acoustic waves at f^^ 6 " 58.8 Hz. 

C . The Cavity Resonance 

A large number of the experimental upward spikes in the experi- 
mental curve in Figure 1 are identified with the theoretical resonance 
frequencies of the standing plane acoustic wave reflected from the 
wooden back panel. A complete list of these cavity resonance frequencies 
is given. 

D. The Acoustic Resonance 

The rest of the experimental upward spikes in the experimental 
curve in Figure 1 are identified with the theoretical resonance fre- 
quencies of the acoustic resonance modes of the panel vibrations. A 
complete list of these acoustic resonance frequencies is given. 

E. The Plate Resonance 

The experimental downward spikes (with the exception of one) in 
the experimental curve in Figure 1 are identified with the theoretical 
resonance frequencies of the plate resonance modes of the panel vib- 
rations. A complete list of these acoustic resonance frequencies is 
given . 

F. The Wooden Back Panel Resonance 

The two most prominent downward spikes in the experimental curve 
in Figure 1 are identified with the theoretical resonance frequencies 
of the "plate” resonance frequencies of the wooden back panel. The first 
sharp downward spike has not been identified before* 
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The detailed features of the experimental curve in Figure 1 have 
been explained by using the six major theoretical aspects of the 
experimental results discussed above, and the theoretical results agree 
with experiment. 
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